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Los  Angeles,  California  90007 

Period:  1  Dec  79  to  30  Nov  80 

1.  Objectives 

The  objective  of  this  study  is  to  explore  and  under¬ 
stand  the  various  conditions  under  which  an  optical  beam 
can  be  generated  that  is  the  phase  conjugate  to  a  given 
input  beam,  when  the  generation  is  mediated  by  nonlinear 
refraction  in  optical  materials.  The  studies  will  explore, 
but  not  be  limited  to,  phase-conjugation  in  liquids, 
glasses,  crystals  and  vapors  in  which  the  optical  beams 
may  be  pulsed  or  cw,  may  be  of  one  or  several  wavelengths, 
and  may  be  free  or  guided. 

2.  Accomplishments  and  Progress 

Accomplishments  and  progress  during  this  reporting 
period  included  the  following: 
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2.1.  We  completed  a  detailed  study  of  optical-beam 
phase- conjugation  by  a  technique  employing  the  photorefrac- 

tive  effect  in  crystalline  BaTiO^.  This  study  also 
included  the  development  of  the  first  theoretical  model 
to  treat  transient  photorefractivity  and  charge  screening 
effects.  Two-beam  effects  and  photoconductivity  measurements 
were  also  used  to  verify  the  theory.  The  results  have 
been  published  as  Publication  3.8  which  is  attached 
as  Appendix  I. 

2.2.  We  have  achieved  the  first  experimental  demon¬ 
stration  of  a  true  time-reversed  vector  replica  of  a  com¬ 
plex  nonuniformly  polarized,  image-bearing,  optical  beam. 

This  is  an  example  of  a  special  kind  of  "phase-conjugation" 
that  cannot  be  achieved  by  holography  or  by  stimulated 
scattering,  but  is  possible  with  four-wave-mixing.  This 
has  been  reported  in  Publication  3.9.  which  is  attached  as 
Appendix  II. 

2.3.  We  have  demonstrated  edge-enhancement  of  an 
optical  image  in  real  time  using  two  different  variations 
of  our  previously  demonstrated  phase-conjugation  technique 
that  employs  the  photo- refractive  effect.  The  results  are 
described  in  Publication  3.10  which  is  attached  as 
Appendix  III. 
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2.4.  Experiments  on  four-wave  mixing  in  photoref ractive 
materials  which  produced  nonlinear  effects  other  than  phase- 
conjugation  were  performed.  The  results  have  been  outlined 

in  Publication  3.11,  a  paper  delivered  at  the  XI  International 
Quantum  Electronics  Conference.  An  abstract  is  attached  as 
Appendix  IV. 

2.5.  We  achieved  the  first  experimental  demonstration 

of  a  wide-angle  (~45°)  narrowband  (~ 0.1cm  optical  filter 

using  phase-conjugation.  Preliminary  results  were  given  at 
XI  International  Quantum  Electronics  Conference,  the 
published  Abstract  for  which  is  listed  as  Pub.  3.5.,  and 
attached  as  Appendix  V. 

2.6.  We  began  a  study  of  the  nonlinear  optical  proce¬ 
dures  of  resonant  iodine  vapor,  tuning  an  argon  laser  through 
several  electronic-vibrational-rotational  resonances  and 
analyzing  the  scattered  radiation  with  a  special  spectro¬ 
meter.  The  first  new  findings  of  this  investigation  will 

be  reported  at  the  December  DEAP-APS  meeting,  and  are 
summarized  in  Publication  3.13  which  is  attached  as 
Appendix  VI. 

2.7.  We  are  continuing  studies  of  resonant  nonlinear 
refraction  and  phase-conjugation  in  sodium  vapor.  Our  results 
include  an  heretofore  unreported  effect,  which  we  call  a 
"tensor  grating".  This  will  appear  first  as  an  Abstract 
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(Publication  3.14)  for  a  forthcoming  talk  at  the  December 
DEAP-APS  meeting.  See  Appendix  VII. 

2.8.  We  have  obtained  the  first  experimental  demon¬ 
stration  of  phase-conjugation  in  the  nanosecond  regime  by 
the  photorefractive  effect  in  BaTiO.^.  Since  this  effect 
"remembers"  the  last  image  conjugated  between  nanosecond 
pulses,  it  has  potential  computer  applications  not  shared 
by  ordinary  four-wave-mixing  effects.  First  results  will 
appear  as  the  Abstract  of  Publication  3.15,  which  is 
attached  as  Appendix  VIII,  and  will  be  given  at  the  Fall, 
1980  meeting  of  OSA. 

2.9.  We  have  achieved  the  first  experimental  demon¬ 
stration  of  a  continuous  wave  phase-conjugating  mirror 
with  gain  (i.e.,  reflectivity  greater  than  unity).  We 
have  used  this  mirror  to  form  self-oscillating  resonators 
and  showed  that  these  resonators  automatically  correct 
phase-distortions  introduced  into  the  resonator.  Details 
of  these  results  given  in  publication  (No.  3.16.),  a 
reprint  of  which  is  attached  as  Appendix  IX.  A  five  minute 
16mm  motion  picture  was  made  of  some  remarkable  self¬ 
oscillation  properties  of  this  mirror  when  faced  with  any 
object  with  a  glint,  such  as  a  kitchen  spatula.  The  movie 
also  shows  the  single-mode  output  from  a  self-oscillating 
cavity  with  aberrator. 
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3.  Publications  from  the  project 

3.1.  "Theory  of  phase-conjugation  in  waveguides  by  four-wave 
mixing,"  R.W.  Hellwarth,  IEEE  Journal  Quant.  Elect.  QE15 , 

101  Feb.  1979. 
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Lett.  33,  404,  Sept.  1978. 
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3 .  Publications  (cont.) 

3.9.  "Generation  of  a  time-reversed  replica  of  a  nonuniformly 

polarized  image-bearing  optical  beam,"  G.  Martin,  L.K.  Lam, 
and  R.W.  Hellwarth,  Optics  Lett.  vol.  5,  pp.  185-187,  May 
1980.  Attached  as  Appendix  II. 
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Jack  Feinberg,  Optics  Letters,  vol.  5,  pp.  330-332,  Aug. 

1980.  Attached  as  Appendix  III. 

3.11.  "Four-wave  mixing  in  photorefractive  materials".  Jack 
Feinberg  and  R.W.  Hellwarth,  Paper  and  Abstract  E.3,  Pro¬ 
ceedings  of  XI  International  Conference  on  Quant.  Elect., 
Boston,  June  23-26,  1980.  Attached  as  Appendix  IV. 

3.12.  "A  wide-angle  narrowband  optical  filter  using  phase-conju¬ 
gation  by  four-wave  mixing  in  a  waveguide",  L.K.  Lam  and 
R.W.  Hellwarth,  Abstract  E.10,  Proc.  XI  Inti.  Quant.  Elect. 
Conference,  Boston,  June  23-26,  1980.  Attached  as  Appendix 
V. 

3.13.  "High  resolution  resonance  Raman  spectroscopy  of  Iodine 
Vapor,"  D.  Kirillov,  and  R.W.  Hellwarth,  Bull.  Am.  Phys. 

Soc.,  DEAP  Meeting,  1-3  Dec.  1980.  Attached  as  Appendix  VI. 

3.14.  "Pulsed  phase  conjugation  due  to  a  tensor  refractive  index 
grating  in  sodium  vapor",  S.N.  Jabr,  L.K.  Lam  and  R.W. 
Hellwarth,  Bull.  Am.  Phys.  Soc.,  DEAP  Meeting  1-3  Dec.  1980. 
Attached  as  Appendix  VII. 
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3.  Publications  (cont.) 

3.15.  "Phase  conjugation  with  nanosecond  laser*  pulses  in  BaTiO^/" 
L.K.  Lam,  T.Y.  Chang,  Jack  Feinberg,  and  R.W.  Hellwarth, 
Abstract,  Bull.  Optical  Soc.,  Fall  1980  meeting,  Chicago. 
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The  following  professional  personnel  were  associated 
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Quantum  Electronics  Conference,  Boston,  June  23. 

5.2.  "Wide-angle  narrowband  optical  filter  using  four- 
wave  mixing  in  a  waveguide",  L.  Lam  and  R.W.  Hellwarth;  paper 
E.10.  XI  International  Quantum  Electronics  Conference, 

Boston,  June  23. 

5.3.  Annual  meeting  of  the  Optical  Society  of  America, 
14-17  October  (Chicago,  Ill.).  "Phase  conjugation  with  nano¬ 
second  pulses  in  BaTiO^",  by  L.  Lam,  T.Y.  Chang,  J.  Feinberg 
and  R.W.  Hellwarth;  paper  Thl3,  16  Oct.,  1980.  Abstract 
published  in  1980  Annual  Meeting  Program  of  OSA. 

5.4.  "Pulsed  phase  conjugation  due  to  a  tensor  refractive 
index  grating  in  sodium  vapor",  S.N.  Jabr,  L.K.  Lam,  and  R.W. 
Hellwarth;  paper  DB8,  Dec.  2,  1980,  American  Phys.  Soc.  -  Div. 
of  Electron  and  Atomic  Physics  Meeting  (at.  U.S.C.). 

5.5.  "High  resolution  resonance  Raman  Spectroscopy  of 
iodine  vapor",  D.  Kirillov  and  R.W.  Hellwarth;  paper  FBll, 
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5.  Talks  presented  (cont.) 


Dec.  2,  1980,  American  Phys.  Soc.  -  Div.  of  Electron  and 
Atomic  Physics  Meeting  (at  U.S.C.) . 

5.6.  Invited:  "Optical  Phase  Conjugation"  by  R.W. 
Hellwarth,  Paper  Q-l,  Dec.  16,  1980.  Lasers  '80  Conference, 

New  Orleans,  La. 

5.7.  "Phase-conjugation  mirror  with  cw  gain",  J.  Feinberg 
and  R.W.  Hellwarth;  paper  Q-4,  Dec.  16,  1980.  Lasers  '80 
Conference,  New  Orleans,  La. 
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Photorefractive  effects  and  light-induced  charge  migration  in  barium  titanate 

Jack  Feinbarg,  0.  Hetman,  A.  R.  Tanguay,  Jr.,  and  R.  W.  Heitwarth 

Electronics  Sciences  Laboratory.  University  of  Southern  California.  Los  Angelos.  California  90007 

(Received  17  September  1979;  accepted  for  publication  1  November  1979) 

W>  propose  a  new  theoretical  model  for  the  light-induced  migration  of  charges  which  mediates 
the  “photorefractive  effect”  (light-induced  refractive  index  change)  in  barium  titanate  and  other 
crystals.  We  also  present  experimental  results  of  various  effects  of  this  light-induced  charge 
migration  in  a  single-domain  crystal  of  barium  titanate,  specifically,  (1)  energy  transfer  between 
two  intersecting  optical  beams,  (2)  optical  four-wave  mixing  and  optical-beam  phase  conjugation, 

(3)  erasure  of  spatial  patterns  of  photorefractive  index  variations,  and  (4)  photoconductivity.  The 
theoretical  model  predicts  the  observed  dependences  of  these  effects  on  (1)  beam  intensities, 
directions,  and  polarizations,  (2)  crystal  orientation,  and  (3)  on  an  externally  applied  dc  electric 
field.  Time  dependences  of  transients  as  well  as  steady-state  magnitudes  are  predicted.  In  this 
model,  identical  charges  migrate  by  hopping  between  adjacent  sites,  with  a  hopping  rate 
proportional  to  the  total  light  intensity  at  the  starting  site.  The  net  hopping  rate  varies  with  the 
local  electric  potential  that  is  calculated  sdf-consistently  from  the  charge  migration  pattern.  In 
barium  titanate  the  charges  are  positive  with  a  density  of  (1.9  ±  0.2)  X 10“  cm*1  at  514  tun.  The 
origin  of  the  charges  and  sites  is  at  present  unknown.  The  hopping  rate  constant  determined  from 
optical  beam  interactions  is  used  to  predict  the  observed  photoconductivity  of  1.3  X 10*'° 
cm  12  *'  W*'  at  514  nm. 

PACS  numbers:  42.65.  -  k,  78.20Jq,  72.40.  +  w,  42.30. Va 


L  INTRODUCTION 

When  light  is  transmitted  through  certain  noncentro- 
symmetric  crystals,  it  causes  a  change  in  the  refractive  index 
which  persists  for  hours  or  longer  in  the  dark  and  can  be 
erased  by  flooding  the  crystal  uniformly  with  light  This 
“photorefractive”  effect  arises  from  a  light-induced  migra¬ 
tion  ar.d  separation  of  charge  in  the  crystal  which  gives  rise 
to  internal  static  electric  fields.  These  fields  produce  refrac¬ 
tive-index  changes  via  the  linear  electro  optic  (Pockels) 
effect. 

We  propose  a  new  theoretical  model  for  the  migration 
of  charges  mediating  the  photorefractive  effect  in  barium 
titanate.  We  also  present  detailed  experimental  studies  of 
charge  pattern  erasure,  and  of  two-wave  and  four-wave  mix¬ 
ing  of  optical  beams  (in  the  milliwatt  range)  in  barium  tita¬ 
nate.  Using  our  theoretical  model  we  are  able  to  predict  the 
observed  dependence  of  wave  mixing  on  the  intensities  and 
polarization  of  the  waves,  and  on  the  wave  directions  relative 
to  each  other  and  to  the  crystal  optical  axis.  Both  the  tran¬ 
sient  and  steady-state  cases  are  discussed. 

It  should  be  pointed  out  that  the  electric  fields  caused 
by  light-induced  charge  migration  are  easily  observable  in 
ferroelectrics  due  to  their  large  electro  optic  coefficients,  but 
such  electric  fields  may  also  be  produced  in  materials  that 
have  inversion  symmetry  and  consequently  lack  a  linear 
electro  optic  effect.  Although  the  fields  would  not  produce 
any  first-order  index  changes  (Pockels  effect)  in  such  materi¬ 
als,  the  fields  could  produce  second-order  index  changes 
(Kerr  effect)  or  energy  level  shifts  (Stark  effect). 

All  of  our  experiments  were  performed  on  a  single-do¬ 
main  2.2  x  2.8  X  4.2-mm  crystal  of  barium  titanate  (Ba- 
TiOj). 1  Between  5  and  1 33  *C,  barium  titanate  is  a  ferroelec¬ 
tric  with  tetragonal  symmetry  C^'.  Our  sample  was  slightly 


wedged  with  faces  cut  approximately  parallel  to  the  (001), 
(010),  and  (100)  planes.  The  crystal  has  a  pale  yellow  color 
and  gives  10*1  extinction  between  crossed  polarizers. 

In  our  “four-wave  mixing”  experiments,  two  beams  of 
light  of  the  same  frequency,  called  “writing”  beams,  inter¬ 
sect  in  the  crystal  and  create  a  periodic  modulation  of  its 
refractive  index.  This  “index  grating”  is  monitored  by  dif¬ 
fraction  of  a  third  “reading”  beam  incident  on  the  grating  (at 
an  angle  that  satisfies  the  Brag*;  condition)  to  produce  a 
fourth  “output”  beam.  If  the  two  writing  beams  and  the 
reading  beam  all  have  the  same  frequency  (and  polarizations 
if  the  sample  is  birefringent),  then  the  Bragg  condition  is 
satisfied  by  a  reading  beam  which  propagates  counter  to  ei¬ 
ther  of  the  writing  beams.  This  geometry  is  identical  to  that 
commonly  used  in  phase  conjugation  experiments,2  and  our 
output  beam  is  observed  to  be  the  phase  conjugate  of  one  of 
the  writing  beams.  However,  in  contrast  to  phase  conjuga¬ 
tion  by  ordinary  nonlinear  refraction,  the  fraction  of  the 
reading  beam  that  is  diffracted  in  the  steady  state  depends 
only  on  the  relative  intensity  of  the  writing  beams  and  is 
independent  of  their  absolute  int  "sity.  A  certain  optical  en¬ 
ergy  (rather  than  power)  must  be  deposited  in  the  crystal  in 
order  to  “write”  a  grating  with  a  given  diffraction  efficiency. 
The  optical  absorption  of  our  sample  was  too  small  to  mea¬ 
sure  (  <  5%),  but  we  estimate  that  several  microjoules  per 
image  element  are  required  to  write  a  grating  of  high 
efficiency. 

If  the  writing  beams  are  turned  off,  the  grating  will  per¬ 
sist  overnight  in  a  darkened  room,  but  reading  the  grating 
erases  it  with  a  decay  rate  which  increases  with  the  intensity 
of  the  reading  beam.  The  grating  can  also  be  erased  by  flood¬ 
ing  the  crystal  with  light  incident  from  an  arbitrary  direc¬ 
tion.  Wavelengths  from  477  to  —900  nm  were  used  to  erase, 
with  the  longer  wavelengths  erasing  at  a  considerably  slower 
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rate  with  the  same  incident  intensity.  A  red  “erasing”  beam 
will  erase  a  grating  previously  formed  by  green  writing 
beams,  and  vice  verse.  Erasing  the  crystal  will  completely 
restore  it  to  its  original  state;  we  have  used  tbe  same  sample 
for  over  a  year  with  no  sign  of  damage  or  discoloration. 
When  making  tbe  grating,  the  rise  time  of  tbe  diffraction 
efficiency  will  decrease  if  the  total  writing  intensity  is  in* 
creased.  For  example,  the  exponential  time  constant  varies 
from  a  few  seconds  to  a  few  milliseconds  as  the  total  writing 
intensity  is  varied  from  10°  to  1  W/cm2.  Both  the  writing 
and  the  erasing  rates  increase  dramatically  when  the  cross¬ 
ing  angle  between  the  writing  beams  is  increased  at  small 

angles 

The  diffraction  efficiency  R  depends  on  the  orientation 
of  the  crystal  with  respect  to  tbe  writing  beams  as  follows.  A 
is  largest  when  the  c  axis  of  the  crystal  is  aligned  normal  to 
the  planes  of  the  intensity  grating,  and  least  when  the  c  axis  is 
parallel  to  these  planes.  A  increases  as  the  square  of  the  inter¬ 
action  length  of  the  beams,  and  we  have  measured  values  of 
R  as  high  as  35%  at  5 15  nm  with  an  interaction  length  (crys¬ 
tal  length)  of  4  mm. 

One  of  the  more  interesting  and  easily  observed  effects 
is  that  the  writing  beams  emerge  from  the  crystal  with  a 
different  relative  intensity  than  when  they  entered  it3-4;  ener¬ 
gy  is  dearly  transferred  from  one  beam  to  the  other.  The 
direction  of  energy  transfer  is  determined  by  the  orientation 
of  thee  axis  relative  to  the  writing  beams.  No  energy  transfer 
is  observed  if  the  sample  is  slightly  vibrated  or  when  the 
beams  do  not  intersect  in  the  sample.  The  magnitude  of  ener¬ 
gy  transfer  can  be  easily  made  to  approach  100%  for  equal 
input  intensities:  one  writing  beam  emerges  from  the  crystal 
almost  extinguished,  while  the  other  writing  beam  emerges 
with  twice  its  original  intensity. 

A  dc  electric  field  externally  applied  along  the  c  axis  of 
the  crystal  enhances  both  the  diffraction  efficiency  and  the 
magnitude  of  the  energy-transfer  effect.  These  effects  will  be 
discussed  in  detail  in  the  following  sections. 

In  Sec.  II  our  model  is  developed,  and  is  compared  to 
our  experimental  results  in  Sec.  III.  Section  IV  applies  our 
model  to  the  experimental  results  of  others,  and  discusses  the 
differences  between  our  model  and  previous  models. 

II.  THEORY 

We  describe  here  a  model  for  charge  migration  in  pho- 
torefractive  crystals  which,  with  two  parameters,  predicts 
the  energy  exchanges  among  two  or  four  optical  beams.  This 
exchange,  in  both  transient  and  steady-state  regimes,  is 
shown  to  depend  on  the  intensities,  polarizations,  and  angles 
of  the  beams,  and  on  an  applied  static  electic  field  for  a  given 
orientation  of  the  crystal.  Given  the  charge-migration  pat¬ 
terns  from  our  model,  one  can  derive  the  quasistatic  electric 
fields  inside  the  crystal,  and  from  this  the  refractive  index 
changes  caused  by  the  electro  optic  effect. 

According  to  our  model,  there  are  a  certain  number  of 
charges  which  can  occupy  a  larger  number  of  sites  in  any  of  a 
large  number  of  permutations.  In  darkness,  each  charge 
stays  fixed  at  a  site,  but  when  a  charge  is  exposed  to  an  opti¬ 
cal  intensity  /,  it  tends  to  “hop”  to  an  adjacent  site  with  a 
probability  per  second  that  is  proportional  to  1.  Let  W,  be 


the  probability  that  a  migrant  charge  occupies  the  nth  site 
which  is  at  position  x,.  If  the  optical  intensity  at  this  site  is 
/,,  then  the  model  may  be  expressed  mathematically  by 

dJV 

~~  -  -  [»'./.  exp(i^.M) 

-  WmIm  expdft^J] ,  (1) 

where  the  sum  is  over  neighboring  sites  m. 

The  rate  constant  Dmm  expf^Sd^,)  measures  the  ten¬ 
dency  toward  light-induced  hopping  from  site  m  to  site  n. 
We  have  written  it  in  terms  of  the  parameters  »  Z>M, 

and  of  the  static  potential  difference  between  the  sites. 

Here  0  equals  q/\kBT ).  where  q  is  the  charge,  ka  is  Boltz¬ 
mann’s  constant,  and  2"  is  the  lattice  temperature.  This  form 
makes  explicit  in  Eq.  (1)  that  the  relative  site  occupation 
probability  in  steady  state  under  weak  uniform  illumination 
will  obey  stadstial  mechanics,  Le.,  that  Wm/WH 
-  expO 3+mm).  Specifically,  <f>nm  is  equal  to  -  6m  where 
=»  4  (x„).  and  ^  (x)  is  the  quasistatic  potential  existing  in 
the  crystal  due  to  internal  charge  migration,  externally  ap¬ 
plied  fields  and  intrinsic  chemical  potentials.  Whenever  light 
is  present,  the  gradient  of  ^  will  cause  a  net  drift  by  hopping 
of  charges  in  time,  away  from  a  stationary  background  of 
neutralizing  charge  (which  we  assume  does  not  hop).  In  our 
model  this  drift  is  governed  by  Eq.  (1),  for  hopping  in  all 
directions. 

In  barium  titanate,  the  omission  from  the  hopping  rate 
of  any  dependence  on  the  probability  of  occupation  of  the 
final  site  is  consistent  with  the  experimental  results.  That  is, 
in  our  experiments  the  site  occupation  probabilities  ap¬ 

pear  to  be  much  less  than  unity,  implying  that  most  sits  are 
unoccupied.  It  is  straightforward  to  account  for  final-site 
occupation  by  appending  factors  \-Wm  to  the  rate  to  hop  to 
site  m  in  Eq.  (1),  should  it  prove  appropriate; 

To  analyze  our  present  experimental  results  we  will 
need  to  consider  only  the  case  where/,  varies  from  the  inter¬ 
ference  of  two  “writing”  optical  beams  of  the  same  temporal 
frequency  whose  complex  electric  field  amplitudes  at  posi¬ 
tion  x  are  e,  is ,  exp((k,-x)and  SZEZ  expOTtj-x),  See  Fig.  l.The 
complex  polarization  vectors  e ,  and  e2  are  normalized  by 
e\4x  =  \,  etc.  Therefore,  we  may  write 


FIG.  I .  Two  writing  beam*  with  wave  vector*  k,  and  k,  intersect  at  an  angle 
28  in  the  crystal  to  produce  an  intensity  grating  with  wavevector 
k  «  k,  —  k..  The  direction  of  the  positive  c  axis  is  shown. 
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f,  «  G0  +  R tG  expOTc-x.) ,  (2) 

where  k  -  k,  -  k*  G  m2E,Ejitfu  and  G^ml,  \ E,  |2  (the 
sum  extends  over  the  two  writing  beams  plus  any  other 
beams  that  are  present).  An  important  parameter  is  the  com¬ 
plex  modulation  index  m  defined  as 

mmG/G0.  (3) 

Note  that  0<|m|<l. 

If  m  <1,  it  is  appropriate  to  try  approximate  solutions 
for  and  d,  of  the  form 

W,  »  W0  +  ReJKexpOlt-x,)  (4) 

and 

4m  -  d*.  +  Re*  exp(i1wO ,  (5) 

where  d,,  is  the  potential  at  x„  due  to  any  charges  external  to 
the  crystal,  applied  fields,  or  chemical  potentials.  Note  that, 
according  to  Poisson’s  equation,  the  potential*  and  charge* 
distribution  amplitudes  4  and  W  must  be  related  by 

d~  Wpq/eejt1 ,  (o) 

where  p  is  the  average  density  (number  per  unit  volume)  of 
sites,  is  the  permittivity  of  free  space,  and  the  screening 
dielectric  constant  is,  in  general, 

€  »£«•*,  (7) 

where  e  is  the  static  dielectric  tensor  and  k  wk /k. 

In  the  process  of  substituting  Eqs.  (2)— (6)  into  Eq.  ( 1)  to 
solve  for  the  amplitude  lV(i)of  the  charge-density  wave,  we 
will  make  some  simplifying  assumptions  that  are  appropri¬ 
ate  for  our  experiments  in  barium  titanate.  First,  the  direc¬ 
tion  of  any  uniform  electric  Add  (  —  graddo)  and  k  are  as¬ 
sumed  to  be  parallel  (or  anti parallel)  to  the  optical  c  axis  in  a 
uniaxial  crystal.  In  this  case  it  will  not  significantly  alter  the 
result  to  assume  that  the  sites  are  equally  (rather  than  ran¬ 
domly)  spaced  by  the  (rms)  distance  /  in  the  z  direction.  Fur¬ 
ther  we  take  Ar/<l,asis  suggested  by  our  data.  In  this  case  the 
form  of  the  result  is  unaltered  by  assuming  that  the  hopping 
rate  is  zero  except  between  nearest-neighbor  sites  when 
it  has  the  value  D.  Typically,  / ~  l  O'4  cm,  and  0  corresponds 
to  room  temperature  (—  l/40eV),  so  that  if  the  internal  elec¬ 
tric  fields  are  below  SO  kV/cm,  $04mn  is  much  less  than  unity 
in  Eq.  ( 1 ).  The  fidd-induced  dectro  optic  effects  we  observe 
suggest  that  in  our  experiments  internal  fields  rarely  ap¬ 
proach  this  value,  and  so  it  is  justifiable  to  keep  only  the 
terms  in  Eq.  (1)  that  are  of  the  lowest  (first)  order  in  4mm- 

With  the  foregoing  approximations,  substitution  of 
Eqs.  (2H6)  in  Eq.  (1)  gives  [by  equating  coefficients  of 
exp(/k.x,)] 


~  *  -  r  [(w  +  mXar2  +  iaf)  +  u>\ . 
dt 

(8) 

Here 

rmDGgkir- 

(9) 

is  a  characteristic  hopping  rate. 

(10a) 

is  the  normalized  charge  wave  amplitude, 

amk/k0 

(10b) 

is  the  grating  wave  vector  normalized  by  and 

/■-(«,-^V(VJ  (10c) 

is  the  uniform  electric  field  strength  normalized  by  a  charac¬ 
teristic  field /q.  The  characteristic  wavevector  k,  is  defined 
by 

k%mpWj/"Jt.T,  (11) 

a  key  parameter  which,  once  determined  by  curve  fitting  to 
experiment,  fixes  the  average  density  p  fVa  of  migrating 
charges.  (In  our  barium  titanate  sample  k»  was  found  to  be 
—0.3k,  which  implies  a  charge  density  — 1014  cm'*.)  The 
characteristic  field in  Eq.  (10c)  is 

f0mk^klT/q.  (12) 

The  characteristic  field  is  seen  to  be  of  order  — 2300 V/an  in 
our  case. 

Equation  (8)  is  a  differential  equation  of  afamiliar 
form.  Even  when  both  the  intensity  /  (r  )  and  the  applied  elec¬ 
tric  field /(/ )  vary  with  time,  the  solution  to  Eq.  (8)  may  be 
written  in  terms  of  simple  integrals.  We  will  use  solutions  of 
Eq.  (8)  in  special  cases  to  analyze  our  experiments  bdow.  In 
these  experiments,  the  charge  distribution  [see  Eq.  (4)]  is 
detected  by  scattering  from  it  an  optical  wave  which  exper¬ 
iences  a  spatial  modulation  of  optical  susceptibility  of  the 
form  ReX  expfilt-x)  due  to  the  dectro  optic  effect.  The  mod¬ 
ulated  optical  susceptibility  is  written  by  convention  in 
terms  of  the  third-rank  dectro  optic  tensor  R  and  the  ampli¬ 
tude  4  of  the  static  potential  distribution  defined  in  Eq.  (5)  as 

X-ieo^CRkK,*.  (13) 

The  optical  dielectric  function  e„  is  diagonal  for  Cartesian 
coordinates  coinciding  with  the  principal  axes.  In  a  uniaxial 
crystal  its  zz  component  equals  n\,  and  itsxx  and yy  compo¬ 
nents  equal  n\  ( na  and/t,  are  the  ordinary  and  extraordinary 
refractive  indices).  In  barium  titanate  (which  has  tetragonal 
symmetry  4mm)  the  nonzero  components  of  the  electro  op¬ 
tic  tensor  are*  =»  23,  =»  R^,mri3  =  8,  and 

Rm  =  Rm  “  820  (in  units  of  10'12  m/V).  Here  the  r(> 
are  the  conventional  contracted  forms  of  the  electro  optic 
tensor.  Also,  ng  =  2.488  and  n,  =  2.424  (at  515  nm),‘  and 
the  dc  dielectric  constants  at  room  temperature  are  106  and 
4300  measured  paralld  and  perpendicular  to  the  c  axis. 

In  our  experiments,  k  was  always  along  the  c  axis,  and 
our  experiments  were  not  sensitive  to  motion  of  charges  ex¬ 
cept  along  the  c  direction.  However,  if  k  and  graddo  were  not 
parallel  to  the  c  axis,  then  Eq.  (8)  for  w  [from  which  follow  4 
by  Eq.  (6)  and  X  by  Eq.  (13)]  is  generalized  as  follows.  The 
right-hand  side  of  Eq.  (8)  becomes  the  sum  over  j  of  three 
terms  (J  =*  x,  y,  z),  each  of  the  same  form  as  in  Eq.  (8),  but  in 
each  of  which  the  symbols  are  altered  in  definition  as  fol¬ 
lows:  D  becomes  Dj  ( j  »  x,y,  z);  kg  becomes  a  function  of 
the  direction  of  k  through  its  dependence  on  e  in  Eq.  (7);  /2 
becomes  If  ( j  —  x,  y,  z )  in  case  the  sites  are  not  spaced  iso¬ 
tropically;  ctj  *  kj/kg  ( 'j  —  x,  y,  z);  and/becomes 
fj  38  —  Xj  -graddo/  f0  (which  are  the  normalized  space  com¬ 
ponents  of  the  externally  applied  electric  field).  We  will  not 
pursue  this  case,  nor  the  further  generalizations  of  Eq.  (8)  to 
any  of  the  following  cases:  kl,  >  1;  Dmm  existing  beyond  the 
nearest-neighbor  pairs;  or  randomly  placed  or  multiple 
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classes  of  sites.  Such  generalizations  do  not  seem  to  be  need¬ 
ed  to  treat  barium  titanate.  We  now  apply  Eq.  (13)  to  the 
erasure  transients  and  steady-state- two-  and  four-wave  in¬ 
teractions  studied  in  our  experiments. 

A.  Grating  erasure  dm*  constants 

Suppose  a  charge  grating  is  formed  by  two  writing 
beams  which  produce  a  given  grating  k,  and  then  are  turned 
off.  leaving  a  given  grating  amplitude  u  expressed  by  Eq.  (8). 
A  weak  reading  beam  is  turned  on  so  that  its  Bragg-scattered 
signal  monitors  the  relative  magnitude  of  the  grating  with¬ 
out  erasing  it  appreciably.  The  crystal  is  simultaneously 
flooded  with  a  stronger  uniform  (m  *  0)  illumination  (not 
phase  matched)  which  then  causes  the  grating  amplitude  w 
to  decay  exponentially  by  Eq.  (8).  The  signal  is  expected  to 
display  an  exponential  decay,  «  |u>|J,  with  decay  rate  A, 
given  by 

A.-lTU+a*).  .  (14) 

Note  that  A,  is  independent  of  any  applied  dc  field.  The 
dependence  of  this  rate  on  a  is  experimentally  verified  in  Sec. 
III. 


B.  Steady-state  four-wave  mixing 

We  consider  here  the  special  form  of  four-wave  mixing 
that  is  also  known  as  transient  volume  holography.7  Two 
optical  beams  cause  charge  migration,  which  in  turn  pro¬ 
duces  a  steady-state  periodic  electric  field  and  consequently 
a  modulation  of  the  optical  susceptibility  [see  Eq.  (13)].  A 
third  “reading”  beam  RetyEj  exp(ik,.x  —  ivr )  is  introduced 
into  the  crystal  and  Bragg  scatters  off  the  optical  susceptibil¬ 
ity  variation.  In  other-words,  with  the  susceptibility  this 
beam  creates  a  polarization  density  in  the  medium  that  is 
equal  to  Rel£3X-f,  exp(/(k,  —  k2  +  k3)-x  —  ivf  ].  The  ratio 
R  of  the  power  radiated  by  this  polarization  density  (with 
field  polarization  e4)  to  the  input  power  of  the  “reading” 
beam  is  (for  R<\) 


R  = 


7—  <% 
4rt3c 


(15) 


where  L  is  the  effective  interaction  length  of  the  beams,  n3  is 
the  refractive  index  of  the  reading  beam,  and  c  is  the  velocity 
of  light.  This  ratio  R  is  called  the  grating  scattering  efficien¬ 
cy.  If  e4  =  e3,  if  the  third  “reading”  beam  is  nearly  antiparal¬ 
lel  to  the  two  writing  beams,  and  if  the  latter  are  unaffected 
by  the  scattering,  then  one  can  write  a  more  general  formula 
for  R,  valid  for  long  interaction  lengths,  in  which  R  can 
approach  unity.1  Attenuation  can  be  accounted  for  in  Eq. 

( 1 5)  by  using  complex  k  vectors  in  the  foregoing  analysis. 

As  an  example,  consider  the  optimum  scattering  con¬ 
figuration  in  barium  titanate  where  the  applied  electric  field 
and  the  grating  k  are  parallel  to  the  crystalline  (?)  axis  (as 
above),  and  the  reading  beam  is  an  extraordinary  ray  and  the 
writing  beams  are  ordinary  rays.  Let  then 


R^B'-L' 


4/, A 


a  +  if 


( / 1  +  /,  +  /j)J  I  1  +  qJ  +  ictf 
where  /,  and  /,  are  the  intensities  of  the  writing  beams  and  I3 
is  the  intensity  of  the  reading  beam  inside  the  crystal.  Note 
that  when  I}  is  much  weaker  than  /,  and/,,  the  above  expres¬ 


sion  for  A  depends  only  on  the  ratio /|//1aadis  independent 
of  the  total  intensity  /,  +  /2  of  the  writing  beams.  When  an 
optical  intensity  attenuation  of  y  cm'1  exists,  L  is  approxi¬ 
mately  [  1  —  exp(  —  yl  )]/y  for  nearly  collinear  beams  inter¬ 
acting  over  the  crystal  length  /.  The  coefficient  3  is  very 
nearly,  in  this  example, 

B=saunJrJj/0(eJi)V4c ,  (17) 

which  has  a  value  —2.4cm'1  for  barium  titanate  at  515  nm. 
There  is  also  a  term  in  3  (much  smaller  in  barium  titanate) 
proportional  to  rl3  when  i3  has  a  component  along  x  or  y.  If 
all  four  rays  were  of  ordinary  polarization,  B  would  be  modi¬ 
fied  by  replacing n3,rJ3(iy£f  by  n£rl3  which (atsmall  writing 
angles)  is  predimed  to  be  3  times  smaller,  from  previous  mea¬ 
surements  of  the  r,j  coefficients,  just  as  we  find  in  See.  QI  by 
measuring  relative  R  values. 


C.  Steady-state  two-wave  coupling 

We  have  seen  how  two  intersecting  optical  beams  of 
wave  vectors  k,  and  ltj  create  a  quasistationary  spatial  wave 
of  optical  susceptibility  variation  with  wave  vector 
ka>k,  —  kj  and  complex  amplitude  X  given  by  Eq.  ( 1 3).  We 
discussed  in  Sec.  II B  how  a  third  beam  scatters  from  this 
refractive  index  grating.  Clearly,  however,  this  grating  can 
also  “scatter”  the  two  beams  which  produce  it  The  main 
effect  of  this  scattering  is,  as  we  show  now,  to  transfer  energy 
from  one  beam  to  the  other,  the  direction  of  energy  flow 
being  determined  by  the  direction  of  thee  axis.  Since  we  have 
been  able  to  obtain  an  analytic  solution  of  Eq.  (1)  for  the 
charge  distribution  in  the  case  where  one  writing  beam  is 
much  weaker  than  the  other  (|m|  <  1),  we  will  assume  here 
that  beam  2  is  much  stronger  than  beam  1  and  hence  negligi¬ 
bly  affected  by  the  two-wave  interaction.*  The  k  component 
of  the  nonlinear  optical  polarization  density  PNL  has  ampli¬ 
tude  4X-e2£2.  Since  X  is  proportional  to<?(=*  2EiE\)  via 
Eq.  (3)  in  Eqs.  (8),  (6),  and  (13),  this  component  is  linearly 
related  to  E,  (as  in  the  normal  linear  polarization  density) 
and  alters  the  dispersion  relation  for  k,(o>)  to  produce  a  new 
k,(a)  -i-  <5  k,.  Assuming  that  both  the  real  and  imaginary 
parts  of  the  (complex)  alteration  5  k,  are  much  smaller  than 
the  unperturbed  k„  one  obtains  immediately  from  Max¬ 
well’s  equations  that  (for  birefringent  or  optically  active 
media) 


k,<5k, 


(18) 


The  direction  of  8  k,  is  usually  determined  by  boundary  con¬ 
ditions  and  is  generally  not  far  from  parallel  to  k,.  Using  the 
steady-state  solution  of  Eq.  (8)  with  Eqs.  (3).  (6),  and  (13) 
one  obtains  the  prediction  of  the  hopping  model  for  the  al¬ 
tered  steady-state  complex  propagation  vector  of  the  weaker 


wave; 


MW, 


iw'-f0  /  g+if  \ 

2c 1  v  1  +  a-  +  ia/l 

X(*;  ]•«.-*-)(«  •#:). 


(19) 


The  most  easily  observed  part  of  the  change  in  the  propaga¬ 
tion  vector  is  the  imaginary  part  of  <5  k,  which  produces  expo¬ 
nential  growth  or  decay  of  the  weaker  wave  through 
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expC  —  2Im(<5  k ,•££. )  where  L  is  the  interaction  length.  The 
wave  will  experience  gain  or  lou  depending  on  the  signs  of/0 
(which  has  the  same  sign  as  the  charge  q)  and  of  the  appro¬ 
priate  elements  of  R. 

Comparing  Eq.  (19)  with  Eq.  ( 17)  shows  that,  when  the 
two  intersecting  beams  are  extraordinary  rays,  the  exponen¬ 
tial  power  gain  coefficient  [i.e.,  —2  Im(<5  k,*i)]  has  a  magni¬ 
tude  (when  k  —  k0  and/—  0)  of  28,  where  B  is  the  coeffi¬ 
cient  of  Eq.  (17)  that  governed  the  four-wave  mixing 
efficiency.  In  barium  titanate  2 B  is  predicted  to  be  —  5  cm'* 
at  515  run,  independent  of  the  intensity  of  the  strong  beam, 
as  was  easily  verified.  The  detailed  dependences  of  gain  on 
beam  angles,  polarization,  and  applied  fields  that  is  predict¬ 
ed  by  Eq.  (19)  have  been  verified  experimentally  as  described 
in  Sec.  Ill 

O.  Moving  crystal 

If  the  crystal  is  moving  with  a  uniform  velocity  v,  then 
one  can  see  that  Eq.  (8)  for  the  amplitude  to  of  the  charge 
grating  is  altered  by  adding  a  tenn  ik-vto  to  the  left-hand  side. 
The  equation  is  stdl  immediately  integrable.  It  is  seen  that 
small  velocities (Ttr/k— 10"*  cm  sec'1) quench  grating  for¬ 
mation.  as  was  observed  experimentally.  Fora  /VI,  a  reso¬ 
nance  in  R  is  predicted  when  w-k  —  —Taf. 

E.  Photoconductivity 

If  a  uniform  potential  gradient  (in  excess  of  any  chemi¬ 
cal  potential)  exists  in  the  uniformly  illuminated  crystal,  and 
Ohmic  contacts  allow  a  one-dimensional  current  flow,  one 
can  calculate  from  Eq.  0)  the  photoresistivity  pt.  One  finds 
that  p,  is  inversely  proportional  to  the  long-wavelength 
(a— 0)  decay  rate  JT  of  w.  This  decay  rate  is  given  under 
special  conditions  in  Eq.  (9),  or  more  generally,  as  explained 
after  Eq.  (13).  Since  this  decay  rate  is  a  (photoinduced)  di- 


I,  I,  I,  BEAMS  t,  l,  CRYSTAL  l, 

MISALIGNED  rotated 


FIG.  2.  Sceady-Mjt*  energy  coupling  between  two  beam*  in  a  poied  single 
crystal  of  barium  man  tie.  The  direction  of  the  positive  c  axis  is  indicated.  In 
(a)  two  bams  with  equal  incident  intensities  (230  rnW/cm1)  emerge  with 
different  imensaita.  In  (b)  the  beam  are  misaligned  so  aa  to  not  intersect  in 
the  crystal,  and  emerge  with  their  intensities  unchanged.  In  (c)  the  crystal  is 
either  repoled  or  simply  rotated  1 80*  to  reverse  the  direction  of  the  positive  e 
aais.  The  direction  of  energy  coupling  also  reverses.  The  spots  were  pro- 
duced  by  etpourg  a  strip  of  photographic  paper  to  the  beams  after  they  had 
passed  through  the  crystal.  No  dc  Acid  waa  applied  (except  prior  to  the 
experiment  in  order  to  pole  the  crystal). 
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FIG.  3.  Experimental  setup  showing  writing  beams  with  intensities  /,  and 
reading  beam,  and  erasing  beam.  For  steady-state  energy-coupling  ex¬ 
periments.  the  intensity  of  writing  beam  1  (after  the  sample)  wan  monitored 
aa  writing  beam  2  waa  blocked  and  unblocked.  For  four-wave  mixing  ex¬ 
periments.  reading  beam  scatters  off  the  grating  produced  by  the  writing 
beams  and  is  deflected  by  beam  splitter  BS,  into  a  photomultiplier  (PMT). 
For  erase- rate  experiment*,  the  wriffhfc  beams  are  both  blocked  and  an 
erasing  beam  is  suddenly  turned  on.  The  decay  of  the  grating  is  monitored 
with  the  PMT  and  analyzed  by  computer.  The  intensity  of  theerasrng  beam 
is  controlled  by  a  Pockeis  cell  followed  by  a  linear  polarizer 


electric  relaxation  rate,  it  is  not  surprising  that  the  result  is 
the  familiar  formula 

p,C€0~r-',  (20) 

where  e  is  as  in  Eq.  (7)  with  k  taken  to  be  the  direction  of  the 
potential  gradient.  At  1  W/cmJ,  r  was  observed  to  be  —  5 
sec'1 ,  (see  Fig.  7)  and  Eq.  (20)  predicts/?, — 5  X 10* (l  cm  (for 
€  =  102),  compared  to  the  measured  resistivity  of p, 

—  8  X  10*  H  cm  with  the  above  illumination. 

III.  EXPERIMENTAL  RESULTS 
A.  Steady -state  energy  coupling 

The  predictions  of  Eq.  (19)  were  experimentally  veri¬ 
fied  by  energy-coupling  experiments,  in  which  two  beams  of 
equal  intensity  are  crossed  in  the  crystal  and  one  is  observed 
to  emerge  with  more  intensity  and  the  other  with  less.  (See 
Fig.  2.)  The  predicted  dependence  of  the  magnitude  of  the 
energy  coupling  effect  on  the  grating  wave  vector  k  (through 
a  =  k  /k0)  was  checked  by  measuring  the  steady-state  cou¬ 
pling  as  a  function  of  the  crossing  angle  26  of  the  two  writing 

Feinbw^waf.  I30t 


FIG.  4.  Utenaiiy  of  writing  beams  *s  a  function  of  time.  One  beam  a  mom- 
tend.  «kib  the  ocher  bam  is  turned  on  at  the  dm  indicattd  by  the  tnt 
enow  aad  turned  off  at  tlw  second  arrow.  Each  bwm  initially  had  an  inten¬ 
sity  of  <00  mW/an1  baton  catering  the  cryMaL  Theriot  tane  of  the  gnting 
is  mu  to  bestow  1/10  mo.  For  banns  with  equal  incident  intensity 
(*•*!).  the  magnicedeof  the  energy  coupled  AJ  into  or  out  of  the  other 
beam  is  sota  to  be  equal  aad  opposite  (hen  about  20%).  The  plots  for  tbr 
two  beams are  slightly  (tiapiaenri  for  clarity. 


beams,  since  k  «  2n0<wc*1  sin0,  where  /*„  is  the  index  of  the 
refraction  appropriate  for  the  direction  and  polarization  of 
the  writing  beams.  A  multimode  argon  ion  laser9  supplied 
two  writing  beams  polarized  perpendicular  to  the  plane  of 
incidence.  The  optical  path  difference  of  the  writing  beams 
from  beam  splitter  1  (see  Fig.  3)  to  the  sample  was  kept 
within  0.3  cm,  well  within  the  approximate  5-cm  coherence 
length  of  the  multimode  laser.  The  beams  were  unfocused 
and  about  2  mm  in  diameter  at  the  sample.  The  interaction 
length  of  the  two  beams  was  determined  by  the  2.2  mm 
thickness  of  the  sample  and  was  independent  of  crossing  an* 
gle  for  the  range  of  angles  used. 

The  exit  intensity /,  of  beam  1  was  measured  with  beam 
2  blocked  (no  grating).  Beam  2  was  then  unblocked,  and  the 
altered  intensity  /,  +  d/„  was  observed  to  be  either  greater 
or  smaller  depending  on  the  orientation  of  the  c  axis  of  the 
crystal  relative  to  the  writing  beams.  (See  Figs.  2  and  4.)  The 
fractional  intensity  changed/,//  was  measured  with  ordi¬ 
nary  waves  for  a  range  of  grating  wave  vectors  (beam  cross¬ 
ing  angles).  The  prediction  of  Eq.  (19)  for  the  steady-state 
fractional  gain  for  AIX  </,  is 


yi  »  -2  lm(8k.xL)=» 


01 

c 


a  +  if 
1  4-  a1  +  iaf 


(21) 


Here  nQ  is  the  ordinary  index  of  refraction,  a  =»  k  /k^ 
f  =  E%vviimi/fo,  and  f0  *  k^kg  T  /q.  A  plot  of  this  function 
and  the  experimentally  measured  values  are  shown  in  Fig.  5. 
The  only  unknown  parameter  is  the  wave  vector  fc<>,  which 
was  found  to  give  the  best  least-squares  fit  to  the  measured 
data  when  ka  —  0.31(2)n0w/c  for  X  »  5 15  nm. 

Energy-coupling  efficiencies  measured  with  both  writ¬ 
ing  beams  at  633  and  at  438  nm  peaked  at  kg**  0.3 \(2)n<fli/c 
at  each  wavelength.  This  implies  that  k#l  is  constant  and, 
from  Eq.  (1 1),  that  the  density  of  filled  trap  sites  pW3  in¬ 
creases  approximately  as  the  square  of  the  frequency  of  the 
writing  beams  over  the  limited  frequency  range  explored. 

Physically,  the  peak  efficiency  of  energy  coupling  was 


FIG.  3.  Fraction  of  steady-state  energy  couplad  d///bccwaenrwcbaameia 
undoped  barium  tjtaaaia  ace  function  of  grating  ware  vector  A.  Tbctrian- 
fies  (d  )  are  the  experimental  data  with  no  externally  applied  dc  electric 
tod,  and  the  tohd  line  is  the  leaet^quanx  fit  to  Eq.(Zl)  with  the  character- 
aac  wave  vector  given  by  fc, — QJ  l(2)*a»/c-  The  open  dreim  (»>  are  the 
experimental  data  with  a  voltage  of  1000  V  applied  acmes  tba  crystal 
(d  »14aa  between  ilenrodee)  parallel  to  than  mm.  Solid  line  lathe 
theoretical  St  to  Eq.  (21)  uaiag  the  same  value  of  k,  and  tba  applied  de 
electric  Add  ineide  the  crystal  ae  the  adjustable  parameter. 


always  found  to  occur  at  the  same  outside  crossing  half-an¬ 
gle  (21*)  of  the  writing  beams  for  all  three  wavelengths  tried. 
This  suspicious  coincidence  was  confirmed  by  determining 
ko  through  other  experiments,  as  follows. 

B.  Eras*  rate 

A  uniform  beam  of  light  (m  =»  0)  incident  on  thesample 
at  an  arbitrary  angle  will  erase  any  grating  that  may  have 
been  previously  stored  in  the  sample,  causing  the  diffraction 
efficiency  to  decay  at  a  rate  .4,  =*  2T(1  +  cr)  where 
r  =*  DGJc%l 1  and  G0  is  proportional  to  the  intensity  of  the 
erasing  beam  [see  Eq.  (14)]. 

To  verify  the  dependence  on  k  (through  a)  of  this  rate, 
the  angle  between  the  writing  beams  was  fixed  while  a  weak 


FIG.  6.  Decay  rate  A.,  a*  a  function  of  erasing  beam  intensity  at  various 
grating  wave  vectors  (crossing  angles).  The  grating  k  in  units  of  nv/c  is  (x) 
k  m  0.03,  (D)  k  »  0.22,  (®)  k  »  0.37,  (V)  k  ■»  0.32.  The  decay  me  is  linear 
in  the  erasing  beam  intensity  for  small  rates,  but  begins  to  saturate  for  rates 
greater  than  I  sec'1. 
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extraordinary  reading  beam  (/^  <  10"3/,*,,),  incident  on 
the  sample  at  the  Bragg  angle,  monitored  the  grating.  Leav¬ 
ing  the  reading  beam  on,  both  writing  beams  were  turned  off 
and  the  erasing  beam  turned  oil  Decay  curves  were  accumu¬ 
lated  for  different  erasing  beam  intensities  /,  and  stored  in  an 
on-line  computer  (Hewlett  Packard  9825A  with  6940B  mul¬ 
tiprogrammer).  The  decay  rate  of  each  curve  was  deter¬ 
mined  and  the  slope  of  a  plot  A ,  versus/,  found.  The  experi¬ 
ment  was  then  repeated  for  different  crossing  angles  of  the 
writing  lv*"w,  thereby  changing  k  (see  Fig.  6).  A  plot  of  the 
slopes  asa  function  of  k  2  gavea  good  straight-line  fit  (see  Fig. 
7)  and  yielded  a  value  for 

kg  *■  0.32(2)it^>/c 

at  A  »  515  nm,  in  excellent  agreement  with  the  k^  deter¬ 
mined  by  energy  coupling  experiments.  From  Eq.  (11)  this 
values  of  k a  gives  a  density  of  filled  trap  sites  of 
pW3  »  1.9(2)X  10“  cm*3.  The  origin  of  these  sites  is  un¬ 
known  and  may  be  due  to  traps  at  impurity  or  defect  sites. 

C.  Sign  of  charge  carriers 

There  are  many  ways  to  determine  the  sign  of  the 
charge  carriers  in  photorefractive  materials.10-1 1  The  sign  of 
the  charge  carriers  in  barium  dtanate  was  experimentally 
determined  to  be  positive  by  energy-coupling  experiments, 
as  follows.  From  Eqs.  (18),  (13),  and  (3),  the  direction  of 
energy  coupling  will  depend  only  on  the  sign  of  the  electro¬ 
optic  coefficient,  the  orientation  of  the  crystal  relative  to  the. 
writing  beams,  and  on  the  sign  of  the  charge  carriers. 

In  order  to  avoid  any  confusion  in  determining  the  sign 
of  the  charge  carriers,  a  method  was  used  which  is  indepen¬ 
dent  of  the  definition^/  the  positive  direction  of  the  e  axis. 

The  crystal  was  first  poled  by  painting  silver  electrodes 
on  the  c-axis  faces,  slowly  heating  it  to  about  3  *C  below  the 
Curie  temperature  (Tc  =*  133  *C),  and  then  applying  a  dc 
poling  field  of  3.5  kV/cm  along  the  c  axis.  The  crystal  was 
then  slowly  cooled  to  room  temperature  with  the  poling  field 
still  applied  to  produce  a  single-domain  crystal.  For  conve¬ 
nience.  we  define  the  positive  c-axis  direction  as  pointing 
toward  the  electrode  that  had  been  connected  to  the  negative 
terminal  of  the  applied  dc  voltage.  (The  final  determination 
of  the  sign  of  the  charge  carriers  is  independent  of  this 
choice.) 

Next  the  sign  of  the  electro  optic  coefficient  was  deter¬ 
mined  to  be  positive  relative  to  the  positive  c  axis  by  placing 
the  crystal  between  crossed  polarizers  with  the  c  axis  orient¬ 
ed  at  45*  to  the  direction  of  polarization.  A  Babinet-Soleil 
compensator,  also  placed  between  the  polarizers,  was  adjust¬ 
ed  to  null  out  the  birefringence  of  the  crystal  at  5 1 5  nm.  A  dc 
field  was  applied  in  the  same  direction  as  the  previous  poling 
field,  and  the  compensator  was  adjusted  to  give  a  new  null 
that  corresponded  to  a  decrease  in  the  difference  of  the  ex¬ 
traordinary  and  ordinary  indices  when  the  applied  field  was 
increased.  From  (13), 

*,(£)*/».  -|n3r„£, 

where  n,(£)  is  the  extraordinary  index  with  applied  field  £, 
nt  is  the  extraordinary  index  with  no  applied  field,  etc.  Since 


FIG.  7.  PVxofUwslopwofFtf.  6  (it,  tb*  dacay  rat*  aoraaiizad  uy  d» 
cratini  born  tnteuicy)***  function  ofk2.  The  drd«n  until*  expariawtul- 
ly  measured  diu  uxl  th*  unifht  line  is  the  beat  Iwst-squires  St.  Tb*  ruin 
of  th*  intercept  to  th*  slop*  of  the  above  graph  yield*  a  vale*  of 
ka  —  0J2(2>».»/c.  The  ui  urccpt  fives*  values  ZT /G„  —  IOmc*1  W*‘  cm * 
used  to  cilnilu*  th*  photoconductivity  in  Eq.  (20). 


|«2  rJS|  >  I*o#jiI,  *  decrease  in  ».(£)  —  «„(£)  with  increas¬ 
ing  E  implies  that  r33  >  0. 

The  poled  crystal  was  then  placed  at  the  intersection  of 
two  optical  beams,  and  the  beam  that  entered  the  crystal  at 
an  acute  angle  to  the  positive  c  axis  was  always  observed  to 
emerge  with  the  weaker  intensity  (see  Fig.  2).  From  Eq.  (19) 
[and  Eq.  (13)  with  r33  >0],  this  implies  that/0 and  hence q  is 
positive. 

In  similar  experiments  on  LiNbOj,  the  direction  of  en¬ 
ergy  coupling  was  opposite  to  that  observed  here.1 1  Assum¬ 
ing  that  the  pyroelectric  determination  of  the  positive  direc¬ 
tion  of  the  c  axis  used  in  that  work12  agrees  with  the 
definition  used  here;  the  signs  of  the  charge  carriers  must  be 
opposite  in  the  two  materials.  However  this  point  should  be 
checked. 

D.  Phase  and  absolute  magnitude  of  index  grating 

The  magnitude  of  the  steady-state  energy  coupling  be¬ 
tween  two  bams  depends  on  the  component  of  the  index 
grating  that  is  90*  out  of  phase  with  the  intensity  distribution 
of  the  two  beams  (i.e.f  on  the  imaginary  part  of  6k,).  Howev¬ 
er,  the  diffraction  efficiency  R  with  which  the  grating  scat¬ 
ters  a  third  reading  beam  depends  on  the  magnitude  of  the 
grating  and  not  on  its  phase  with  respect  to  the  writing 
beams.  Therefore,  by  separately  measuring  the  fraction  F  of 
two-beam  energy  coupling  and  the  diffraction  efficiency  de¬ 
scribed  in  Sec.  II  B,  real  and  imagina-y  parts  of  Sk  in  Eq.  (19) 
can  be  determined.  From  Eq.  (19) 'with  no  applied  electric 
field  (/=  0),  the  real  part  of  Sk  should  be  zero,  and 

F=2R  1,2 

when  all  beams  have  the  same  polarizations,  m  =  1,  and 
R  <  1.  The  above  relation  was  experimentally  verified  within 
10%  error  over  a  range  of  crossing  angles  and  confirmed 
that  the  index  grating  with  wave  vector  k  is  90*  ±  10*  out  of 
phase  with  the  incident  intensity  distribution. 

Equation  (19)  for  the  fraction  of  energy  coupling 
reaches  a  peak  at  k  =  k0and  predicts  1 3%  coupling  for  a  2.2- 
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APPLIED  VOLTAGE  (KV) 

FlO.  &-  Four-wave  muting  diffraction  efficiency  Rttt  function  of  applied 
voltage  at  fixed  grating  spacing  (*a«0.0S  Wc).  The  volute  m  applied 
peralM  to  the  c  uu*.  The  crystal  wn  2.3  nun  in  width  between  the  e-uis 
faces.  The  aroaaw  ate  theexparimenul  data  and  the  solid  line  iatha  theoreti¬ 
cal  fit  to  Eq.(  16)  (if  the  internal  applied  field  wen  40%  of  the  applied 
voltage  divided  by  the  crystal  width). 


mm  crystal  length  using  ordinary  polarizations. 11  This  com¬ 
puted  value  is  36%  less  {Jhan  the  peak  measured  value  of  0.28 
with  mat  1  in  Fig.  5,  which  may  be  due  to  contributions  from 
higher-order  terms  in  Eq.  (8)  when  m  is  near  unity. 


E.  Energy  coupling  and  diffraction  efficiency  with  an 
applied  electric  field  ( f  #0) 


To  confirm  the  predictions  of  our  theory  in  the  presence 
of  an  applied  field,  both  the  diffraction  efficiency  R  and  the 
amount  of  energy  coupled  between  the  two  writing  beams 
were  measured  with  and  without  +1000  V  applied  across 
the  crystal  parallel  to  the  +•  c  direction,  and  as  a  function  of 
the  crossing  angle  of  the  writing  beams. 

From  Eq.  (8>  the  photoinduced  electric  field  with  peri¬ 
odicity  k  is 


£(,)  -  R.  ^ |  ”**■»>)  • 

(22) 

The  diffraction  efficiency  R  is  proportional  to  |£  j1  [see  Eq. 
(16)]  and  has  the  form 


R  <x 


0*  +  /’ 


(23) 


(l+oV+O2/2’ 

The  dependence  of  R  on  applied  voltage  is  shown  at  a 
fixed  crossing  angle  in  Fig.  8.  This  curve  represents  a  vertical 
slice  through  the  curves  shown  in  Fig.  9  and  shows  that  R  is 
approximately  quadratic  in /for  (i.e.,  a<l).  The  fit  to 

the  measured  data  is  excellent 

If  there  were  an  intrinsic  field  in  the  BaTiO,  crystal,  as 
suggested  by  Chen,'3  then  the  slope  of  the  graph  would  be 
nonzero  when  the  applied  field  was  zero.  The  above  data  and 
the  observed  fact  that  R— *0  as  a— -0  puts  an  upper  limit  on 
such  a  field  of  <  300  V/ctn.  However,  some  assym- 

metry  was  observed  in  R  for  electric  fields  applied  parallel 


and  antiparallel  to  thee  axis  and  indicates  that  some  intrinsic 
field  may  be  present  in  our  sample. 

The  results  of  energy  coupling  experiments  are  shown 
in  Fig.  3.  There  is  a  noticeable  shift  in  the  peak  of  the  curve 
toward  smaller  crossing  angles  as  the  applied  electric  field  is 

increased,  and  an  increase  in  the  magnitude  of  energy  coup¬ 
led,  especially  at  small  crossing  angles.  The  best  fit  to  the 
deta  gave  a  value  of  the  applied  field  which  was  73%  of  the 
actual  field  applied.  This  discrepancy  is  possibly  due  to  the 
poor  contact  made  by  the  silver  paint  electrodes  or  by 
boundary  layers.14 

IV.  DISCUSSION 

We  here  discuss  the  main  differences  between  our  the¬ 
ory  and  previous  theories  and  experiments  on  barium  rita- 
nate  and  similar  photo-refractive  materials.  A  semiperman¬ 
ent  light-induced  change  was  first  observed  in  lithium 
niobate  and  in  other  ferrodectrics  by  Ashkin.13  It  soon  be¬ 
came  apparent  that  this  photorefractive  effect  or  “optical 
damage”  was  caused  by  a  redistribution  of  charges  in  the 
material  creating  an  electric  field,  which  then  altered  the 
index  of  refraction  through  the  electro  optic  effect.13  It  was 
first  postulated  theta  large  intrinsic  electric  field  was  neces¬ 
sary  to  enable  thecharges  to  drift  in  the  material,'3  but  it  was 
subsequently  demonstrated  that  diffusion  alone  could  ac¬ 
count  for  the  observed  effects.16  High-quality  holographic 
storage  was  subsequently  demonstrated  in  these  materials 
using  two  optical  beams  of  one  wavelength  to  create  a  vol¬ 
ume  hologram,  and  a  beam  of  a  different  wavelength  to  read 
it17  Many  theories  have  been  proposed  to  describe  various 
aspects  of  volume  holograms  in  photorefiractive 
materials. I*-23 

An  important  goal  of  these  theories  is  to  predict,  from  a 
given  incident  intensity  distribution  /(x),  the  resulting  re¬ 
fractive  index  modulation  n(x).  The  approach  followed  in 
previous  papers  was  to  first  calcualte  the  number  of  charge 


FIG.  9.  Theoretical  plots  from  Eq.  (13)  of  the  four-weve  mixing  direction 
efficiency  R  as  a  function  of  grating  wave  vector  k  for  various  tppiied  elec¬ 
tric  fields.  The  numbers  on  each  plot  are  the  values  of /(the  applied  dc 
electric  field  in  units  of/,  *  T/q). 
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carriers  excited  to  the  conduction  band,  and  then  to  compute 
their  motion  from  the  forces  of  drift  and  diffusion.  In  all  of 
these  theories,  the  number  of  charge  carriers  excited  into  the 
conduction  band  was  always  taken  to  be  proportional  to 
/(x).  This  assumption  ignores  the  depletion  of  charge  carri¬ 
ers  available  to  be  excited  in  the  high-intensity  regions  of  the 
sample  and  is  not  valid  except  when  writing  with  short  opti¬ 
cal  pulses.  Our  theory  shows  that  the  density  of  filled  sites  is 
signiAcantly  altered  by  charge  migration  as  steady  state  is 
approached,  even  for  very  small  modulation  index  m,  and 
cannot  be  assumed  constant  (For  a  »  t,  and  m  »  1,  the 
modulation  of  the  filled  site  density  approaches  unity.)  Con¬ 
sequently.  in  our  theory,  we  include  bopping  rates  that  are 
proportional  to  the  product  of  J(x)  and  the  probability  that 
the  site  at  x  is  filled. 

In  addition  most  of  the  previous  theories  On  charge 
transport  limit  themselves  to  the  initial  stage  of  grating  for¬ 
mation  and  avoid  the  effects  of  the  photoinduced  field  itself 
in  the  equation  of  motion  for  the  charges.24  However,  it  is 
just  this  field  created  by  charge  bopping  that  counters  the 
force  of  diffusion  and  enables  the  system  to  reach  equilibri¬ 
um.  This  is  described  explicitly  in  Eq.  (1)  and  results  in  the 
steady  state  as  well  transient  solutions  that  are  evident  from 
Eq.  (8)  and  is  an  important  characteristic  of  our  theory. 

A  key  parameter  in  our  theory  is  the  wave  vector  k^  [see 
Eq.  (11)].  The  value  of  k 0  will  vary  with  material  and  from 
sample  to  sample  according  to  the  density  of  filled  trap  sites 
available  for  hopping  when  excited  by  light  of  a  given  fre¬ 
quency.  Alphonse  et  aJ.13  have  noted  a  diffraction  efficiency 
in  lithium  niobate  which,  as  a  function  of  grating  wave  vec¬ 
tor  k,  has  a  distinct  turning  point  similar  to  our  Fig.  3  and 
which  they  note  is  “a  striking  experimental  departure  from 
(their)  theory."  Applying  our  theory,  the  peak  of  their  plot 
yields  a  value  of  k0  =  0.36(7)n<u/c  at  488  nm  and  a  filled  trap 
density  of  7x  I0‘s  cm'3  at  the  optical  wavelength  employed. 
This  is  quite  close  to  the  value  reported  here  for  our  sample 
of  undoped  barium  titanate.  Data  reported  by  Townsend 
and  LaMacchia'9  on  barium  titanate  also  show  R~k2  for 
small  k  and  a  roll  off  at  about  k0  *  0.2 nca/c  at  488  nm  which 
implies  a  filled  trap  density  of  7  X  IQ15  cm'3.  Two  data  points 
of  R  versus  k  reported  by  Peltier  and  Micheron  in  a  sample  of 
bismuth  silicon  oxide23  imply  a  filled  trap  density  somewhat 
greter  than  2X  10“  cm'3.  If  the  trap  sites  are  due  to  impuri¬ 
ties  or  defect  sites,  it  is  not  clear  why  the  filled  trap  densities 
should  be  similar  within  a  factor  of  2  in  two  different  samples 
of  barium  titanate  and  in  a  sample  of  iron-doped  lithium 
niobate.  Further  experiments  may  perhaps  provide  some 
physical  insight  as  to  why  the  concentration  of  filled  carriers 
is  so  similar  in  different  materials,  and  why  the  density  of 
traps  appears  to  vary  as  the  square  of  the  optical  frequency  of 
the  writing  light. 

In  summary,  we  have  presented  a  new  charge-transport 
model  to  describe  the  formation  of  volume-index  gratings  in 
photorefractive  materials.  Predictions  of  the  model  for  the 


diffraction  efficiency  of  the  grating  and  for  the  fraction  of  the 
energy  coupled  between  the  two  writing  both  in  the 
transient  and  steady  state,  are  well  confirmed  by  experiment. 
A  closed-form  solution  is  presented  for  values  of  the 

optical  modulation  m  (and  which  fits  the  dots  well  for  rrtsx  1) 
and  predicts  the  maximum  in  scattering  efficeincy  R  versus 
grating  wave  vector  k  observed  by  us  and  by  others.  This 
model  also  predicts  the  observed  charge- pattern  erasure  rate 
as  a  function  of  k  and  the  observed  dependence  of  diffraction 
efficiency  on  the  intensities  and  polarizations  of  the  optical 
beams.  The  model  also  predicts  the  observed  photoconduc¬ 
tivity  and  predicts  effects  of  sample  velocity. 
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We  demonstrate  the  generation  of  the  time-reversed  replica  of  an  incident  monochromatic  image-bearing  optical 
beam  that  ia  nonuniformly  polarised.  A  form  of  degenerate  four-wave  mixing  in  liquid  CSj  of  beam*  at  532  nm  is 
employed.  We  discuss  source*  of  distortion  in  the  vector  replica  using  standard  theory  of  nonlinear  beam  interac¬ 
tion. 


Introduction 

The  generation  of  a  time-reversed  replica  (Le.,  the  phase 
conjugate)  of  a  monochromatic,  image-bearing  optical 
wave  has  been  demonstrated  previously  by  a  variety  of 
techniques,  but  only  for  cases  in  which  the  wave  has  a 
known  and  uniform  (usually  linear)  state  of  polariza¬ 
tion.1  Here  we  report  the  generation  of  a  complete 
time-reversed  vector  replica  (or  phase  conjugate)  of  an 
incident  monochromatic  image-bearing  optical  wave 
whose  state  of  polarization  has  been  altered  by  a 
nonuniformly  birefringent  window  so  that  it  is  randomly 
nonuniform  over  the  wave  front. 

We  employ  a  kind  of  degenerate  four- wave  mixing 
process  in  which  the  incident  wave  is  made  to  overlap 
two  counterpropagating  plane  pump  waves  whose 
electric  Held  vectors  are  rotating  (circularly  polarized) 
in  opposite  directions.  Theory  predicts  that  the  wave 
generated  by  the  resulting  third-order  nonlinear  po¬ 
larization  in  an  isotropic  medium  is  the  desired  time- 
reversed  vector  replica,  provided  that  the  incident  and 
pump  waves  are  mutually  temporally  coherent  where 
they  overlap  in  the  nonlinear  medium  (CSj).2 

The  earliest  demonstration  of  the  generation  of  a 
time-reversed  replica  of  a  uniformly  polarized  optical 
wave  (i.e.,  scalar  phase  conjugation)  was  by  Kogelnik, 
who  used  standard  holographic  methods  of  still  pho¬ 
tography.3  The  first  demonstration  of  nearly  instan¬ 
taneous  scalar  phase  conjugation,  by  Stepanov  et  al., 
employed  a  degenerate  four-wave  mixing  process  that 
was  described  as  a  kind  of  transient  holography.4 

Our  experimental  arrangement  for  conjugating  a 
nonuniformly  polarized  wave  is  described  below.  This 
arrangement  is  similar  to  one  proposed  by  Zel’dovich 
and  Shkunov5  for  conjugating  a  wave  front  of  uniform 
but  arbitrary  polarization.  The  theory  of  this  vector- 
image  conjugation  is  outlined  in  the  section  on  theory, 
and  the  area  of  applicability  of  this  process  is  described 
in  the  discussion  section. 


Experiment 

The  experimental  arrangement  we  have  employed  to 
study  vector-wave  phase  conjugation  is  shown  sche¬ 
matically  in  Fig.  1.  A  linearly  polarized  optical  pulse 
of  ~I5-nscec  duration  and  ~1-MW  peak  power  at 
532-nm  wavelength  is  generated  by  a  frequency-doubled 
Nd:  YAG  laser.  This  wave  is  split  into  three  beams  by 
beam  splitters  BSl  (35%  reflecting)  and  BS2  (50^  re¬ 
flecting)  of  Fig.  1.  The  two  beams  from  BS2  form 
counterpropagating  pump  waves  G  and  H,  which  are 
circularly  polarized  (counterrotating)  by  quarter-wave 
plates  Ql  and  Q2  and  focused  into  the  interaction  region 
at  C  by  lenses  Ll  and  L2.  The  100%  reflecting  dielectric 
mirrors  Ml  and  M2  do  not  affect  the  polarizations  of  G 
and  H. 

The  wave  that  is  split  off  by  beam  splitter  BSl  is 
linearly  polarized  at  angle  <t>  to  the  vertical  by  polarizer 
P2  and  expanded  by  telescope  BE  to  a  diameter  of 
about  2  cm.  This  beam  F  illuminates  the  object  0  and 


Fig.  1.  Schematic  of  apparatus  for  observing  (at  plane  S) 
backscattered  images  of  incident  beam  F  from  region  C.  The 
backscattering  is  caused  either  by  an  ordinary  mirror  at  C  or 
by  nonlinear  mixing  of  the  incident  image-bearing  beam  F 
with  counterpropagating  plane  pump  beams  G  and  H  in  liquid 
CS2  (contained  in  a  2-mm-Iong  cell  at  C).  The  role  of  the 
various  optical  elements  is  discussed  in  the  text. 
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is  than  directed  by  mirror  BS3  through  phase  distorter 
Dl  and  polarization  distorter  D2  and  is  focused  by  lens 
L3  into  region  C  in  the  nonlinear  medium  (CS?),  where 
it  crosses  the  counterpropagating  pump  beams  at  an 
angle  of  about  4°.  All  three  beams  (F,  G,  and  H)  have 
a  diameter  of  —0.2  mm  in  the  cell,  where  they  overlap 
for  the  entire  2-mm  length  of  the  cell. 

The  phase-conjugate  image  is  observed  behind  beam 
splitter  BS3  and  is  recorded  on  photographic  print 
paper  at  observation  plane  S,  which  is  approximately 
the  same  distance  from  the  interaction  region  as  is 
object  0- 

Figure  2  compares  photographic  records  at  plane  S 
of  waves  backscattered  from  the  incident  image-bearing 

I  I 


Fir.  2.  Backscattered  images  recorded  at  plane  S  (Fig.  1). 
Toe  images  of  column  I  were  backscattered  from  beam  F  by 
an  •.-rdinary  plane  mirror  at  C  (the  focal  plane  of  lens  L3  of  Fig. 
1 1.  The  images  of  column  II  were  backscattered  from  beam 
F  by  nonlinear  mixing  with  the  counterpropagating  plane 
pursp  waves  in  liquid  CS?  (as  depicted  in  Fig.  1).  Various 
cor.dit  ions  of  phase  and  polarization  distortion,  as  well  as  in¬ 
cident  and  backscattered  beam  polarizations,  result  in  the 
irr.aies  shown  here  in  rows  a-e,  as  is  explained  in  the  text. 
The  backscattered  fraction  from  nonlinear  mixing  was  of  the 
order  of  1*7.  These  images  were  recorded  on  FSC  polycon- 
tra.-:  G  photographic  printing  paper.  Images  in  column  I  are 
averages  of  a  few  shots.  Images  in  column  II  are  averages  of 
ah'  ut  in  shots. 


beam  F  either  by  a  plane  mirror  at  position  C  (column 
I)  or  by  the  phase-conjugation  process  when  the  cell  of 
CSg  is  at  C  (column  II).  The  pictures  of  rows  a,  b,  and 
c  of  Fig.  2 -were  taken  without  polarizer  P3  under  the 
following  conditions: 

Picture  la  is  the  mirror-reflected  image  of  object  0 
.(a  bent  wire)  with  the  distorting  plates  Dl  and  D2  re¬ 
moved.  (The  image  is  inverted  because  of  the  presence 
of  focusing  lens  L3.)  Picture  Ila  is  the  phase-conjugate 
image  with  the  CSo  ceil  at  position  C  and  distorters  Dl 
and  D2  absent.  Pictures  lb  and  lib  are  the  corre¬ 
sponding  images  but  with  phase  distorter  Dl  in  place. 
Pictures  Ic  and  lie  are  the  corresponding  images  when 
both  phase  distorter  Dl  and  polarization  distorter  D2 
are  in  place.  The  cross  superposed  on  this  image  arises 
from  polarization  distorter  D2,  which  consisted  of  four 
quadrants  of  quarter- wave  retarding  plastic,  each  of  a 
different  orientation.  Pictures  Id,  le,  lid,  and  lie  are 
taken  with  linear  polarizer  P3  placed  as  indicated  in  Fig. 
1:  first,  in  row  d,  with  its  axis  aligned  parallel  to  the 
initial  polarization  set  by  P2  of  image-bearing  beam  F, 
and  then,  in  row  e,  aligned  90°  to  this  direction. 

Comparison  of  the  pictures  in  column  I  of  Fig.  2  with 
those  in  column  II  shows  that  our  arrangement  for  phase 
*  conjugation  not  only  restores  phase  distortions  created 
in  an  image-bearing  beam  but  also  corrects  polarization 
scrambling.  To  check  the  degree  of  the  polarization 
restoration  of  image-bearing  beam  F,  we  oriented  the 
F-beam  polarizer  P2  at  different  angles  4>  to  the  vertical 
and  in  each  case  measured  the  fraction  of  the  phase- 
conjugate  wave  polarized  orthogonal  to  the  input  di¬ 
rection  (e.g„  we  measured  the  ratio  of  intensities  in 
pictures  lie  and  lid).  We  found  the  ratio  to  be  less  than 
1:100  for  all  cases.  Beam  restoration  was  observed  to 
deteriorate  as  the  pump -beam  polarizations  were  al¬ 
tered  from  circular. 


Theory 

Our  experimental  conditions  are  expected  from  theory 
to  produce  a  quite  accurate  time-reversed  vector  replica 
of  a  nonuniformly  polarized  incident  wave  of  amplitude 
Fj(x),  i  =*  x,yj.  This  can  be  seen  from  the  expression 
in  esu  (valid  for  small  replication  efficiencies)  for  the 
vector  amplitude  E,(x),  i  =*  x,y,z,  of  the  backscattered 
wave  at  positions  x  in  front  of  the  region  of  nonlinear 
beam  interaction2: 

£,(x)  =  12t  iw(nc)-'LT,jCjktmFi,*(x)GtHm.  (I) 

Here  the  summation  over  repeated  space  indices  0, 

/,  m)  is  assumed.  The  waves  are  ail  at  angular  fre¬ 
quency  w  in  a  nonlinear  medium  having  refractive  index 
n  and  nonlinear  susceptibility  tensor  c;*/m  evaluated  at 
the  four  frequency  arguments  (-w,  -w,  w,  w).  The 
incident  wave  £*  (x)  overlaps  two  counterpropagating 
pump  waves  (whose  vector  amplitudes  G  and  H  in  the 
interaction  region  have  spatial  components  G,  and  Hi) 
over  an  effective  interaction  length  L.  The  tensor  T„ 
yields  the  projection  of  any  vector  onto  the  plane 
transverse  to  the  F  and  E  beams. 

For  isotropic  media  and  optical  wavelengths  such  ns 
we  employ,  the  c,,ki  tensor  has  the  form 

Cijkt  *  abijbkt  +  b6,k0ji  +  c <$, ify*. 


(2) 
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where  the  coefficients  a,  b,  and  c  may  exhibit  temporal 
and  spatial  dispersion.  When  the  nonlinear  polariza¬ 
tion  arises  from  a  nonpropagating  physical  change,  such 
as  the  electronic  and  molecular  reorientation  mecha¬ 
nisms  that  will  predominate  in  our  CSj  medium,  then 
symmetry  requires  that  b  »  c.  For  CSj,  a  variety  of 
existing  measurements  gives  values  (appropriate  to  532 
nm)  4»c~4X  10“13  esu  and  a  —  8  x  10“u  esu  to 
within  10*0,  provided  that  electrostrictive  and  thermal 
effects  can  be  neglected.6 

If  the  angle  9  between  beams  F  and  G  (see  Fig.  1)  is 
much  less  than  unity,  and  if  | b  —  c |  «  |6  +  c|,  then  one 
finds  from  Eqs.  (1)  and  (2)  that  having  G  and  H  oppo¬ 
sitely  circularly  polarized  gives  the  desired  result  that 
E,  <*  F,  *  to  order  02.  Of  course,  a  real  beam  will  not  be 
perfectly  circularly  polarized,  so  assume  for  the  purpose 
of  evaluating  Eq.  (1)  that 

G  m  (£r  +  dL)G  (3) 

and 

H  -  (iL  +  84R)H,  (4) 

where  4r  and  4l  are  the  complex  unit  vectors  repre¬ 
senting  right  and  left  circular  rotations,  respectively,  of 
the  optical  electric  vector  [equaling  (i'  ±  i$')/2l/ 2  for 
waves  directed  along  ±2']  and  the  fractions  t  and  b  of 
“wrong”  polarization  are  assumed  to  be  much  less  than 
unity.  Substituting  Eqs.  (2M4)  in  Eq.  (1)  gives  (in 
vector  notation)  for  the  baclcscattered  field 

E  *  12ffiw(nc)-1Z,[(a  +  +  %c)F*  +  Ya/K] GH. 

(5) 

Here  the  error  field  K  is  defined  as  £Fy  -  yFx  so  it  is 
orthogonal  to  the  desired  conjugate  part  F* .  The  error 
coefficient  f  is 

/  *  (6  +  c)  |(sin2  0  +  (1  +  cos2  9)(e  +  5)jF,*Fy* 

+  i(o  - 1)  cos  9  (F,*F,*  -  Fy*Fy*)\J 
r-  -  (b  -  c)  cos  d  (FX*FX*  +  Fy*Fy*)/F 2,  (6) 

in  which  we  have  omitted  terms  higher  than  first  order 
in  t  and  6.  In  Eq.  (5),  terms  in  the  coefficient  of  F*  that 
are  fractionally  smaller  by  terms  of  the  order  of  «,  5, 02, 
and  (5  —  c)/b  have  also  been  neglected. 

The  power  fraction  |//(6  +  c  +  2a)|2  of  the  “wrong” 
polarization  of  the  backscattered  beam  is  seen  to  have 
terms  of  the  order  of  (1)  the  power  fractions  t2  and  52  of 
wrong  polarization  in  the  pump  beams,  (2)  the  broken 
symmetry  ratio  |(6  —  c)/(b  +  c)|2,  and  (3)  sin4  9,  as  well 
as  the  three  cross  terms  of  the  order  of  the  square  root 
of  the  product  of  any  two  of  the  above.  It  is  seen  to  be 
experimentally  feasible  to  keep  the  fraction  of  wrong 
polarization  generated  (at  any  place  x  in  the  phase 
front)  well  below  1%,  as  we  observed. 

Discussion 

One  can  interpret  the  three  terms  in  the  nonlinear 
susceptibility  Eq.  (2)  acting  in  Eq.  (1)  to  produce  the 
backscattered  field  E  as  a  superposition  of  four  gTatings 
of  linear  refractive  index  off  which  the  beams  G  and  H 
scatter  to  generate  E.  The  four  gratings  are  (1)  beam 
G  scattering  off  the  scalar  grating  F*  •  H  to  give  a  term 
E  <*  cGF*  •  H,  (2)  beam  H  scattering  off  the  scalar 


grating  F*  •  G  to  give  ^term  E  «  b  HF*  •  G,  (3)  beam  G 
scattering  off  the  tensor  grating  F*H,  and  (4)  beam  H 
scattering  off  the  tensor  grating  F*G,  the  last  two  giving 
a  term  E  «  a  F*H  •  G.  Consider  grating  (1)  of  ampli¬ 
tude  «  F*  •  H  and  of  wave  vector  —2 nus/c.  The  grating 
must  not  move  by  an  appreciable  fraction  of  the  grating 
period  ir c/rtu  during  the  longest  characteristic  time  te 
associated  with  the  pulse  shape  and  the  physical 
mechanism(s)  responsible  for  the  grating  (mainly  mo¬ 
lecular  reorientation  in  CSz,  for  which  tr  —  I  psec). 
This  is  achieved  if  the  frequencies  of  the  F  and  H  beams 
are  the  same  to  within  ~tc-1  (or  if  they  are  correlated 
in  a  special  way). 

The  grating  also  must  not  move  by  an  appreciable 
fraction  of  its  period  during  the  time  it  takes  the  reading 
beam  G  to  propagate  over  the  length  L  of  the  interaction 
region.  This  is  achieved  if  the  frequencies  of  the  F  and 
H  beams  are  the  same  to  within  — c IL.  In  either  case, 
the  frequency  v  of  the  reading  beam  G  can  be  quite 
different  (hundreds  of  wave  numbers)  from  the  fre¬ 
quency  of  F  and  H  without  spoiling  the  phase  matching 
for  the  finite  lengths  L  of  practical  interest.  However, 
the  backscattered-image  shrinkage  factor  (w/v)  may 
become  a  problem.  All  these  considerations  arise  in  the 
usual  scalar-wave  phase-conjugation  experiments, 
whether  with  photographic  or  transient  grating  media, 
and  apply  to  the  other  three  gratings  me-  'ioned 
above. 

Aside  from  the  obvious  application  to  restoring  a 
beam  reflected  back  through  a  birefringent  distorting 
medium,  vector-wave  replication  may  be  used  in 
imaging  microscopic  inclusions  in  such  a  medium. 
Inclusions  that  have  a  complex  (nonreal)  optical  di¬ 
electric  tensor  do  not  sustain  time-reversed  propagation, 
and  so  they  may  be  imaged  by  a  lens  at  the  observation 
plane  (S  of  Fig.  1)  free  of  the  distortion  caused  by  the 
surrounding  inhomogeneous  medium,  which  has  a  real 
(though  anisotropic)  dielectric  tensor. 

In  conclusion,  we  have  demonstrated,  as  was  expected 
from  theory,  the  generation  by  a  form  of  four- wave 
mixing  of  a  time-reversed  replica  of  a  monochromatic 
image-bearing  optical  wave  of  arbitrary  nonuniform 
polarization. 
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Edge  enhancement,  a  type  of  optical  image  processing,  is  performed  in  a  photorefractive  materiel  in  real  time  and 
with  low  incident-light  intensities  (10~J  W/c m2).  We  calculate  the  expected  images  using  two  different  four-wave 
mixing  geometries,  which  show  good  agreement  with  the  images  that  we  experimentally  observe  using  a  single-do¬ 
main  crystal  of  BeTiOs  aa  the  photorefractive  material. 


In  this  Letter  we  exploit  the  nonlinearity  inherent  in 
the  photorefractive  effect  to  demonstrate  and  explain 
edge  enhancement  of  an  optical  image.  Edge  en¬ 
hancement,1  aa  well  aa  other  kinds  of  optical  image 
processing,2  has  been  performed  previously  in  photo¬ 
refractive  materials.  These  materials  (BaTi03, 
LiXb03,  BSO,  and  BGO,  to  list  a  few)  respond  to  low- 
light  levels  (—1  mJ/cm2)  to  create  an  efficient  volume 
hologram  that  will  persist  for  many  hours  in  the  dark.3 
F or  our  experiments  we  use  a  poled  crystal  of  undoped 
barium  titanate  (BaTiCh),  which  measures  2.2  mm  X 
2.S  mm  X  4.2  mm  and  is  used  with  no  externally  applied 
electric  field  or  electrical  contacts. 

Two  different  geometries  are  used  to  produce  two 
different  kinds  of  edge  enhancement.  In  the  first,  a  real 
image  of  the  object  is  focused  in  the  sample,  and  in  the 
second  the  spatial  Fourier  transform  of  the  object  is 
imaged  in  the  sample.  In  either  case,  the  object-bearing 
beam,  or  object  beam,  intersects  a  reference  beam  in  the 
sample  to  make  an  intensity  interference  pattern, 
which,  by  the  photorefractive  effect,  creates  a  semi¬ 
permanent  index  grating  (or  volume  hologram)  in  the 
crystal.  A  third  beam,  called  the  reading  beam,  inci¬ 
dent  at  the  Bragg  angle  of  this  index  grating,  is  deflected 
by  the  grating  and  is  separated  for  observation  by  a 
beam  splitter  (see  Fig.  1)  to  form  a  real  image  of  the 
original  object.  This  image  is  a  faithful  replica  of  the 
original  object  only  if  the  intensity  of  the  object  beam 
is  less  than  the  intensity  of  the  reference  beam.  If  this 
condition  is  violated,  then  the  image  will  show  markedly 
enhanced  edges  (see  Fig.  2). 

To  understand  the  origin  of  edge  enhancement,  let 
the  optical  electric  fields  of  the  incident  beams  be  of  the 
form 

E,(x)  -  Re(£;(x)#;  exp‘ik>’«  “  '-*>], 

where  £/(x),  4j,  and  k;  are  the  electric  field  amplitude, 
the  polarization  vector,  and  the  wave  vector,  respec¬ 
tively.  of  the  two  writing  beams  0  ■  1,2)  and  of  the 
reading  beam  O'  »  3),  all  measured  in  the  crystal  and  all 
with  frequency  w.  Consider  the  special  case  in  which 
all  three  incident  beams  are  plane  waves  so  that  the 
£;(x»  are  uniform  across  the  beams.4  The  uniform 
writing  oeams  will  create,  where  they  intersect,  a  peri¬ 
odic  array  of  index  grating  planes  in  the  crystal  with 
wave  vector  k  »  -  kj.  If  the  crystal  is  aligned  such 

that  its  optic  axis  (c  axis)  is  parallel  to  the  grating  wave 


vector  k,  then,  from  the  theory  of  the  photorefractive 
effect,5  the  steady-state  index  grating  will  have  an  in- 
tensity-scattering  efficiency  R  given  by 

R  »  |££Cnt|2  (1) 

for  a  reading  beam  incident  at  the  Bragg  angle  and  for 
R  «  1.  Here  L  is  the  interaction  length.  The  coeffi¬ 
cient  B  is  2.4  cm-1  for  barium  titanate  at  515  nm  with 
an  extraordinary  reading  beam.  C  is  a  geometrical 
factor  of  the  order  of  unity.  The  key  factor  in  Eq.  (1) 
is  the  modulation  depth  m,  which  is  given  by 

m  a  2EiE2*6\  •  e2*/(]£i|2  +  IE2I2  +  IE3I2).  (2) 

In  this  experiment,  the  intensity  of  the  reading  beam 


Fig.  1.  Experimental  setup  showing  writing  beams  with  in¬ 
tensities  1 1  and  1 2  (ordinary  polarisation)  and  reading  beam 
with  intensity  (extraordinary  polarization).  The  angle  0 
measured  outside  the  BaTi03  crystal  is  18°.  The  direction 
of  the  c  axis  of  the  crystal  is  as  shown.  Adjustable  neutral- 
density  filters  (F)  are  used  to  vary  the  intensities  of  the  three 
incident  beams.  M  are  mirrors,  and  BSi  and  BS2  are  50% 
reflecting  beam  splitters.  BSj  is  a  4%  reflecting  beam  sp'itter. 
The  lens  (L3)  and  the  pinhole  (H)  spatially  filter  the  object 
beam,  which  is  recollimated  by  lens  L2  and  focused  in  the 
crystal  by  lens  L».  The  polarizer  (P)  passes  the  signal  beam 
but  prevents  scattered  light  from  the  writing  beams  from 
reaching  the  camera  (Cl  or  the  optical  multichannel  array  (not 
shown).  The  two  unfocused  optical  beams  measure  ap¬ 
proximately  2  mm  in  diameter  at  the  crystal,  which  has  a 
thickness  l  -  2.2  mm. 
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Fig.  2.  Images  of  a  moustache  comb  showing  the  second  type 
(Fourier)  of  edge  enhancement  described  in  the  text  In  the 
photograph  on  the  left,  the  intensity  of  the  object  beam  is 
1/100  the  intensity  of  the  reference  beam.  In  the  photograph 
on  the  right,  the  intensity  of  the  object  beam  is  100  times  the 
intensity  of  the  reference  beam.  The  teeth  of  the  comb  are 
separated  by  0.8  mm.  The  intensity  doublets  described  in  the 
text  and  shown  in  Fig.  4  are  barely  resolvable  here  at  the  edges 
of  the  teeth  of  the  edge-enhanced  comb. 

is  made  much  less  than  the  intensity  of  either  of  the 
writing  beams,  the  writing  beams  are  linearly  polarized 
parallel  to  each  other  and  are  ordinary  rays  in  the 
crystal,  and  the  reading  beam  is  linearly  polarized  or¬ 
thogonally  to  the  writing  beams  and  is  an  extraordinary 
ray  in  the  crystaL  This  choice  of  polarizations  for  the 
three  incident  beams  both  maximizes  the  ratio  of  signal 
to  scattered  noise  and  minimizes  the  interaction  be¬ 
tween  the  reading  beam  and  the  writing  beams  during 
grating  formation.  Under  these  conditions,  Eq.  (2)  for 
m  reduces  to 

m  -  2(AA)1/2/(A  +  A),  (3) 

where  /i  and  A  are  the  intensities  of  the  object  and 
reference  beams,  respectively.  Note  that  m  is  now  in¬ 
dependent  of  the  total  intensity  of  the  object  and  the 
reference  beams  (A  +  A),  and  therefore  the  reflectivity 
R  will  depend  only  on  their  relative  intensity  (A/A)- 
The  type  of  edge  enhancement  depends  on  how  the 
object  beam  is  focused  in  the  crystal,  as  follows.  In  the 
first  method  of  edge  enhancement,  an  object  (such  as 
a  binary  transparency  or  a  resolution  chart)  is  placed 
in  the  object  beam,  and  a  lens  (Lx  in  Fig.  1)  is  adjusted 
to  produce  an  inverted,  demagnified,  real  image  of  the 
object  in  the  crystal.  Let  this  image  have  a  (now  non- 
uniform)  intensity  A(x).  If  A(x)  varies  slowly  com¬ 
pared  with  the  grating  spacing  2ir/|k|,  then  m  can  be 
taken  as  a  local  modulation  index  m(x), 

m(x)  *  2(71(x)A)1/2/(AU)  +  A)-  (4) 

If  the  object  beam  is  more  intense  than  the  reference 
beam,  then  for  bright  regions  of  the  object,  A(x)  »  A, 
and  mix)  will  be  much  less  than  its  maximum  value  of 
unity.  Similarly,  for  dark  regions  of  the  object,  I i(x) 
«  A,  and  mix)  will  still  be  small.  However,  at  some 
point  xo  in  the  transition  region  between  a  light  and  a 
dark  region  in  the  object  (i.e.,  at  an  edge),  A(*o)  *  A, 
which  gives  m(xo)  equal  to  unity,  and  consequently  in 
the  region  around  xo  a  local  grating  with  a  large  dif¬ 
fraction  efficiency  is  produced.  The  reading  beam  will 
be  preferentially  scattered  by  this  grating,  and  the 
corresponding  image  will  therefore  have  that  edge  en¬ 
hanced.  The  very  bright  and  the  very  dark  regions  of 
the  object  are  relatively  ineffective  in  producing  a 
grating;  only  at  an  edge  will  the  intensity  of  the  object 


beam  match  the  intensity  of  the  reference  beam  to 
produce  a  grating  with  a  large  diffraction  efficiency. 

Figure  3  is  a  plot  of  the  calculated  and  experimentally 
observed  images  when  the  object  is  a  single  slit  oriented 
perpendicularly  to  the  plane  of  Fig.  I.  The  theoretical 
curves  were  generated  by  assuming  that  the  intensity 
in  the  crystal  is 

A(x)  »  A  rect(x/A/d),  (5) 

where  rect(X)  *  1  if  |X|  <1  and  rect(X)  »  0  otherwise, 
d  is  the  width  of  the  slit,  and  M  is  the  magnification  of 
the  focusing  lens  (Li  in  Fig.  1).  The  above  expression 
for  th&object-beam  intensity  is  substituted  into  Eq.  (4), 
the  resulting  expression  for  mix)  is  substituted  into  Eq. 
(1),  and  a  plot  of  R  versus  x  is  calculated.  Note  the 
enhancement  of  the  edges  of  the  image  as  the  ratio  I \JI 2 
is  increased  in  Fig.  3. 

If  instead  of  a  binary  transparency  the  object  has 
shades  of  gray,  then  that  shade  of  gray  will  be  enhanced 
that,  when  illuminated  by  a  strong  object  beam,  has  an 
intensity  that  matches  the  intensity  of  the  reference 
beam.  An  intensity  contour  map  can  be  generated  in 
this  way  by  simply  varying  the  intensity  of  the  reference 
beam. 

In  the  second  method  of  edge  enhancement,  instead 
of  forming  a  real  image  of  the  object  in  the  crystal,  the 
spatial  Fourier  transform  of  the  object  is  imaged  in  the 
sample  plane.  Consider  an  object  composed  of  a  single 
slit  In  the  crystal,  the  electric-fieid  distribution  of  the 
object  beam  is  the  spatial  Fourier  transform  of  the  slit 
(i.e.,  the  familiar  single-slit  diffraction  pattern).  If  A 
»  1 2,  then  in  the  crystal  the  central  maximum  0f  the 
Fourier  transform  will  have  an  intensity  greater  than 
the  intensity  of  the  reference  beam  and,  by  Eqs.  (1  >  and 
(5),  will  produce  a  localised  grating  with  a  poor  dif¬ 
fraction  efficiency.  Toward  the  wings  of  the  diffraction 
pattern,  the  intensity  of  the  object  beam  is  comparable 
with  the  intensity  of  the  reference  beam,  and  conse¬ 
quently  a  high-efficiency  grating  will  be  formed.  Re¬ 
construction  of  the  image  with  the  reading  beam  and 
lens  Li  is  equivalent  to  an  inverse  Fourier  transform. 
Notice,  however,  that  the  low  spatial  frequencies  of  the 
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Fig.  3.  Theoretical  and  experimentally  observed  plots  of  the 
intensity  of  the  observed  image  of  a  single  slit  0.8  mm  wide 
when  a  real  image  (demagnification,  1/6)  of  this  object  is  fo¬ 
cused  in  the  sample  by  a  100-mm  focal-length  lens.  The 
theoretical  plots  have  been  smoothed  to  take  into  account  the 
finite  resolution  of  the  optical  detection  system.  The  ratio 
of  the  intensity  of  the  object  beam  to  the  intensity  of  the  ref¬ 
erence  beam  (L/A)  is  (a)  0.5,  (b)  5,  (c)  50. 
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object,  which  are  contained  in  the  central  lobe  of  the 
Fourier  transform,  have  been  attenuated,  while  the  high 
spatial  frequencies  have  been  enhanced. 

Figure  4  is  a  plot  of  the  theoretical  and  experimentally 
observed  intensity  profiles  in  the  image  plane  when  the 
object  is  a  single  slit  3  mm  wide.  A  200-ram  focal-length 
lens  is  used  (Li  in  Fig.  1)  and  positioned  one  focal  length 
away  from  both  the  object  and  the  crystal  to  produce  the 
Fourier  transform  of  the  object  in  the  crystal  The 
theoretical  curves  are  computed  by  taking  the  Fourier 
transform  of  a  single  slit,  multiplying  the  amplitude  of 
each  frequency  component  by  a  weighting  factor  pro¬ 
portional  to  m  in  Eq.  (4),  recombining  the  terms,  and 
then  squaring  them  to  give  the  intensity  profile.  That 
is,  the  object  beam  was  considered  to  have  an  optical 
electric-field  distribution  given  by 

Ei(x)  *  Ei  sin  ( vxfd ),  (6) 

where  d  is  the  width  of  the  object  slit  Then 

I  /*40/d 

/oWd(*)b  I  £i(*0  •£2/(|£l(X/)|*+|E2|2) 

I  J-AO/d 

cos  2vx'xdx'  I  .  (7) 


where  a  spatial  cutoff  frequency  of  40  d-1  is  used  to  take 
into  account  the  finite  width  of  the  interaction  region 
of  the  crystal,  and  the  computed  intensities  are 
smoothed  by  averaging  adjacent  points  to  duplicate  the 
finite  resolution  of  the  detection  system.  In  the  mea¬ 
sured  intensity  profiles,  the  enhancement  of  the  high- 
spatial-frequency  components  (edges)  is  seen  to  increase 
while  the  low-spatial-frequency  components  disappear 
as  the  intensity  of  the  object  beam  is  increased  relative 
to  the  reference  beam.  Note  the  appearance  of  inten¬ 
sity  doublets  at  the  edges,  which  is  characteristic  of  this 
method  of  edge  enhancement. 

By  misaligning  the  object  beam  it  is  possible  to  en¬ 
hance  edges  preferentially  along  one  direction  only. 
For  example,  when  the  object  beam  is  adjusted  to  be 
slightly  out  of  the  plane  (by  ~0.5  mm)  defined  by  the 
beams  in  Fig.  1,  the  horizontal  edges  of  the  image  be¬ 
come  prominent,  since  the  high  spatial  frequencies  from 
these  edges  extend  into  the  plane  of  the  other  incident 
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Fig.  4.  Theoretical  plots  of  the  electric  field  and  intensity 
in  the  image  plane  and  the  experimentally  measured  intensity 
using  the  Fourier  geometry  described  in  the  text.  The  ratio 
of  the  intensity  of  the  object  beam  to  the  intensity  of  the  ref¬ 
erence  beam  ih/Ii)  ia  (a)  0.06,  (b)  6,  (c)  60. 


Fig.  5.  Images  of  a  resolution  chart  showing  enhancement 
of  only  horizontal  edges  as  a  result  of  misaligning  the  object 
beam  vertically  by  about  1  mm.  The  radial  pattern  is  1  cm 
in  diameter.  The  ratio  of  intensities  of  the  object  and  refer¬ 
ence  beams  is  1/1000  in  the  photograph  on  the  left  and  1000/1 
in  the  photograph  on  the  right 

beams  and  overlap  with  them  in  the  crystaL  However, 
almost  all  the  vertical  information  of  the  object  will  miss 
the  interaction  region,  and  the  reconstructed  image  will 
lack  vertical  detail  (see  Fig.  5). 

It  should  be  stressed  that  edge  enhancement  is  pos¬ 
sible  only  because  the  diffraction  efficiency  R  of  the 
crystal  is  not,  in  general  a  linear  function  of  the  incident 
writing-beam  intensities.  If,  however,  the  reading  beam 
has  an  intensity  much  greater  than  the  intensity  of  both 
of  the  writing  beams  and  is  left  on  during  the  process  of 
grating  formation,  then,  from  Eq.  (2),  m  will  always  be 
much  less  than  unity  and  R  will  be  approximately  bi¬ 
linear  in  Jj  and  I2.  It  was  experimentally  verified  that 
edge  enhancement  vanished  when  the  reading  beam  was 
made  too  strong,  even  with  /i »  / 2. 

In  conclusion,  we  have  demonstrated  and  explained 
real-time  edge  enhancement  of  optical  images  using  two 
different  geometries.  The  observed  images  are  in  good 
agreement  with  the  theoretical  image  intensities  cal¬ 
culated  from  the  theory  of  the  photorefractive  effect. 

The  author  wishes  to  thank  D.  Heiman,  R.  W.  Hell- 
warth,  and  T.  Strand  for  helpful  conversations.  This 
research  was  supported  by  the  National  Science 
Foundation  under  grant  no.  ENG78-04774  and  the  U.S. 
Air  Force  Office  of  Scientific  Research  under  grant  no. 
78-3479. 
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FOUR-WAVE  MIXING  IN  PHOTOREFRACTIVE  MATERIALS 


Jack  Feinberg  and  R.W.  Hellwarth 
Electronics  Sciences  Laboratory 
University  o£  Southern  California 
Los  Angeles,  California  90007 
(213)  741-6402 


Abstract 


Extensive  results  on  time-dependent  optical  phase  conju¬ 
gation  and  image  processing  by  degenerate  four-wave  mixing  in 
photorefractive  materials  are  correlated  using  a  new  model 
of  light- induced  charge  hopping. 
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POOR-WAVE  MIXING  IN  PHOTOREFRACTIVE  MATERIALS 
Jack  Feinberg  and  R.W.  Hellwarth 
Electronics  Sciences  Laboratory 
University  of  Southern  California 
Los  Angeles,  California  90007 

SUMMARY 

In  ordinary  holography  the  interference  pattern  of  two 
optical  beams  (a  reference  beam  and  an  object  beam)  is  stored 
on  film  and  later  viewed  with  a  third  optical  beam.  In  order 
to  faithfully  reproduce  the  image  of  an  object,  the  intensity 
of  the  object  beam  must  be  made  less  than  the  intensity  of  the 
reference  beam.  In  transient  volume  holography  using  photore- 
fractive  materials  (e.g.,  BaTi03,  BSO,  LiNb03)  interesting  new 
effects  are  observed  when  the  above  condition  is  violated.  Using 
four-wave  mixing  in  certain  geometries  we  have  studied  transient 
holographic  effects  in  which  we  have 

a.  enhanced  any  edges  present  in  the  original  object 
while  perserving  the  broad  overall  features  of  the  image 

b.  performed  edge  enhancement  instead  by  sever ly 
suppressing  the  low  spatial  frequencies  of  the  object;  any  large 
bright  regions  become  dark  and  only  the  high  spatial  frequency 
edges  remain 

c.  demonstrated  real  time  pattern  recognition  using 
four-wave  mixing 

d.  generated  real  time  contour  maps  of  an  object. 


Four -Wave  Mixing  In . 

Jack  Feinberg  &  R.w.  Hellwarth 
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We  describe  a  new  "charge  hopping"  model  which  explains  the 
above  effects  in  both  transient  and  steady  state.  We  report 
experiments  which  elucidate  the  physical  mechanism  of  the 
light- induced  charge  hopping. 

Specifically,  we  have  investigated 

a.  the  number  density  and  sign  of  the  charges  that 
hop  in  response  to  light 

b.  charge-hopping  rates  in  directions  parallel  and 
perpendicular  to  the  optic  axis 

c.  the  dependence  of  diffraction  efficiency  (of  the 
light-induced  hologram)  on  the  total  energy  delivered  to  the 
material,  and 

d.  a  predicted  resonance  in  charge-grating  formation 
when  the  crystal  is  translated  with  a  certain  velocity. 
Applications  of  four-wave  mixing  in  photorefractive  materials 
to  phase-conjugation  devices  will  be  discussed. 
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Abstract 


A  highly  convergent  broadband  pulse  ~532  nm  is  optically 


filtered  by  four-wave  mixing  in  an  optical  waveguide.  The 
narrow  passband  is  centered  about  the  frequency  of  guided 
counter-propagating  pump  waves. 
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A  Wide-Angle  Narrowband  Optical  Filter  Using  Phase-Conjugation 
By  Four-Wave  Mixing  In  A  Waveguide* 

L.K.  Lam  and  R.W.  Hellwarth 
Electronics  Sciences  Laboratory 
University  of  Southern  California 
Los  Angeles,  California  90007 

Summary 

When  an  "input"  optical  beam  of  frequency  v  is  launched 
into  a  waveguide  in  which  counter-propagating  monochromatic 
"pump"  waves  exist  at  frequency  ca,  then  a  phase-conjugate 
wave  (pew)  is  generated  at  frequency  2oj-v  by  the  nonlinear 
polarization  in  the  guide  medium.  Depending  on  the  pump 
beam  powers,  the  power  of  the  pew  may  be  larger  than  that  of 
the  input  beam,  but  it  is  essentially  independent  of  the 
angle  of  the  input  be—m  with  respect  to  the  guide  axis,  and 
it  falls  off  rapidly  when  |v-wj  becomes  greater  than  the 
inverse  of  the  propagation  time  through  the  guide.  Using 

i 

these  properties,  we  have  constructed  and  demonstrated  a 
narrowband  optical  filter  having  a  wide  acceptance  angle. 

The  method  is  as  follows.  A  Nd:YAG  laser  pulse,  mainly  in 
a  single  mode  and  of  duration  -20  ns  and  power  -5  MW,  is 
passed  through  a  frequency-doubling  crystal  and  than  a  mixing 
crystal  to  generate  beams  of  doubled  and  tripled  frequencies. 

The  355  nm  beam  pumps  a  Coumarin  dye  cell  in  a  weakly  dispersive 
cavity  to  produce  a  wideband  laser  pulse  centered  at  532  nm. 


A  Wide-Angle  Narrowband  Optical 
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This  wideband  dye  laser  pulse  is  beam-expanded  and  focussed 
as  the  input  beam  into  an  optical  waveguide  using  F/l  optics. 
The  waveguide  is  a  glass  tube  of  I.D.  0.2  mm,  -20  cm  long, 
and  filled  with  CS2  liquid.  The  frequency-douoled  Nd:YAG 
pulse  is  split  into  two  beams  and  weakly  focussed  into  each 
end  of  the  optical  waveguide.  The  pew  emerges  counter  to  the 
dye  laser  beam  and  is  separated  with  a  50%  beam  splitter. 

The  spectral  widths  of  the  input  dye  laser  beam  and  of  the 
output  phase- con jugate  beam  are  measured  with  a  Fabry-Perot 
interferometer.  Filter  bandwidths  are  of  the  order  of 
0.1  cm”*.  Such  an  optical  filter  carbines  the  advantages 
of  the  narrow  bandwidth  of  a  Fabry-Perot  interferometer  and 
the  wide  acceptance  angle  of  a  grating  monochromator. 
Comparisons  of  observed  filter  tr-.nsmission  efficiency 
and  bandwidth  with  theory  are  presented. 

*We  gratefully  acknowledge  the  support  of  the  Air 
Force  Office  of  Scientific  Research  under  grant  No.  78-3479 
and  of  the  National  Science  Foundation  under  Grant 


ENG  78-04774. 


Abstract  Submitted 


/^ppeiidi'x  JE 


Physical  Review 
Analytic  Subject  Index 
Number  33.2QF,  33.20K 


for  the _ DEAP _ Meeting  of  the 

American  Physical  Society 
1-3  Dec.,  1980 


Bulletin  Subject  Heading 
in  which  paper  should  be  pis 
Laser  Spectroscopy _ 


High  Resolution  Resonance  Raman  Spectroscopy  of 
Iodine  Vapor.*  D.  KIRILLOV  and  R.W.  HELLWARTH.  Unlv.  of 
So.  Calif.— We  have  used  the  precise  resonances  of  a 
tuned,  single-mode  S145A  argon  laser  with  the  P(13)  and 
R(15)  rotational  lines  of  the  12o_+(v,,*0)-^H^u(v,»43) 
transitions  of  the  iodine  molecule  to  obtain  the  first 
measurements  of  the  anhaxmonicity  splitting  of  the  Q(13) 
and  Q(1S)  Raman  branches.  The  inelastic  scattering  from 
12  vapor  (-200  mT)  was  collected  within  4°  of  the  for¬ 
ward  direction  and  analyzed  by  a  Fabry- Perot  interfero¬ 
meter  ganged  with  a  grating  spectrometer.  Although  our 
instrumental  resolution  of  ~80Hiz  is  the  highest  resolu¬ 
tion  yet  reported  in  spontaneous  Raman  spectroscopy,  our 
recorded  linewidths  for  forward  scattering  were  still 
instrumental .  The  splittings  which  we  observed  in  the 
Q-lines  were  consistent  with  known  anharmonicity  para¬ 
meters.  We  nave  also  first  observed  that,  whereas  the 
same  resonance  Raman  lines  observed  at  90°  scattering 
have  full  Doppler  linewidths  (~50QMHz) ,  the  lines  near 
backward  scattering  were  narrower  than  instrumental  and 
exhibited  hyperfine  structure. 

*We  acknowledge  the  support  of  Air  Force  Office  of  Sci¬ 
entific  Research  under  grant  No.  78-3479  and  of  the  U.S. 
DOE  through  Lawrence  Livermore  Laboratory  Subcontract 
No.  7509103. 
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Pulsed  Phase  Conjugation  Due  to  a  Tensor  Refractive 
Index  Grating  in  Sodium  Vapor.*  S.N.  JABR,  L.K.  LAM  and 
R.W.  HELLWARTH,  Univ.  of  So.  Calif.— We  observe  the 
scattering  of  a  light  beam  at  frequency  w  from  a  refrac¬ 
tive  index  grating  formed  by  beating  two  other  beams, 
also  at  w,  near  the  one-photon  D-line  resonances  in  so¬ 
dium  vapor.  There  are  two  components  of  comparable  in¬ 
tensity  but  different  spectra.  The  first  'scalar'  com¬ 
ponent,  studied  previously  by  Wandzura,*-  and  Humphrey, 
et  al.,2  arises  from  a  grating  of  excited  sodium  atoms. 
Here  we  examine  a  second  tensor*  component  which  arises 
from  oriented  but  not  optically  pimped  ground-state 
atoms,  and  is  present  even  when  the  two  gTating-forming 
beams  are  orthogonally  polarized  (and  hence  can  produce 
no  grating  of  excited  atoms) .  We  observe,  as  we  ex¬ 
pected,  that  the  scattering  by  this  'tensor'  grating 
occurs  mainly  at  the  Dj  line  when  the  grating  is 
formed  by  nitrogen-laser-pumped  dye-laser  pulses.  We 
present  a  sinple  theoretical  account  of  the  observed 
effects. 

*Work  supported  by  the  AFOSR  grant  No.  79-0098  and  Law¬ 
rence  Livermore  Subcontract  No.  7509105. 
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Photorefractive  materials  have  recently  been  used  to  perform  optical 
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Image  processing  on  a  millisecond  time  scale  *  .  We  report  the  generation  of 
phase-conjugate  waves  ("real-time"  holography)  using  a  single,  ten-nanosecond 
laser  pulse  in  a  crystal  of  the  photorefractive  material  barium  titanate  (BaTiO-j). 

We  find  that  the  optical  energy  necessary  to  form  the  hologram  is  roughly 
independent  of  the  power  of  the  incident  optical  beams  Over  nine  orders  of 
magnitude.  The  application  to  high-speed  optical  computing  will  be  discussed. 

1)  J.  White  and  A.  Yariv,  "Real-Time  Image  Processing  via  Four-Wave  Mixing 
in  a  Photorefractive  Material , 11  Appl .  Phys.  Lett.,  July  (1980). 

2)  Jack  Felnberg,  "Real-Time  Edge  Enhancement  Using  the  Photorefractive  Effect," 
Optics  Letters,  August  (1980). 
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We  demonstrate  a  p hue-conjugating  mirror  that  has  a  continuous-wave  power  reflectivity  much  greater  than  unity 
(gain  ~100).  This  mirror  usea  nonresonant  degenerate  four-wave  mixing  in  a  single  crystal  of  barium  titanate 
(BaTiOs).  With  our  mirror  we  have  (1)  observed  cw  self-oscillation  in  an  optical  resonator  formed  by  this  mirror 
and  a  normal  mirror,  (2)  demonstrated  a  cw  oscillator  that,  in  spite  of  phase-distorting  material  placed  inside  the 
resonator,  will  always  smit  a  TEMoo  mode,  and  (31  demonstrated  an  optical  image  amplifier.  This  mirror  will  work 
at  any  visible  wavelength  and  with  weak  (milliwatt  or  weaker)  pump  beams. 


F has* -conjugating  mirrors  wore  demonstrated  previ¬ 
ously  with  reflectivities  that  are  greater  than  unity,  and 
seif-oscillation  observed,  but  only  for  a  few  nanosec¬ 
onds.1  The  largest  reflectivity  reported  to  date  for  a 
continuous-wave  (cw)  phase-conjugating  mirror  is  only 
17% 2  In  those  experiments,  either  resonant  degenerate 
four-wave  mixing  was  necessary,  which  permitted  op¬ 
eration  only  over  a  small  frequency  range  (~1  GHz),  or 
beams  of  megawatt  power  were  needed.  In  this  Letter, 
we  report  the  first  known  demonstration  of  a  cw 
phase-conjugating  mirror  with  reflectivity  greater  than 
unity.  We  employ  degenerate  four-wave  mixing  of 
milliwatt  beams,  mediated  by  the  photorefractive  effect 
in  a  single  crystal  of  barium  titanate  of  2.2-mm  X 
2.8-mm  X  4.2-mm  dimensions  at  room  temperature. 
The  effect  is  nonresonant  and  operates  over  a  large 
fraction  of  the  visible  spectrum.  The  main  disadvan¬ 
tage  of  this  phase  conjugating  (pc)  mirror  is  its  relatively 
slow  response  time,  of  the  order  of  1  sec  at  the  nominal 
milliwatt-power  levels  of  common  lasers.  (However, 
this  response  time  shortens  inversely  with  the  pump- 
beam  power.) 

By  using  our  pc  mirror,  we  have  (1)  observed  cw 
self-oscillation  in  an  optical  resonator  formed  by  this 
mirror  and  a  normal  mirror,  (2)  demonstrated  wave- 
front  correction  when  a  phase-distorting  medium  is 
placed  inside  the  self-oscillating  resonator  (that  is,  the 
pc  mirror  alters  the  transverse-mode  structure  of  the 
resonator  to  compensate  automatically  for  any  phase 
distortions  in  the  cavity),  and  (3)  demonstrated  optical 
image  amplification. 

To  understand  the  operation  of  this  mirror,  consider 
two  optical  beams,  with  wave  vectors  kL  and  k->,  having 
nonorthogonal  polarizations  and  the  same  angular 
frequency  w.  Call  these  beams  the  writing  beams. 
Where  they  intersect  in  the  crystal,  they  form  an  in¬ 
tensity-interference  pattern  with  wave  vector  k  *  kk  ■ 
ki  —  k2.  Electrical  charges  (of  unknown  origin)  migrate 
in  the  crystal  from  the  peaks  into  the  troughs  of  the 
intensity-interference  pattern  and  eventually  reach  a 
static-charge  distribution.  These  charges  create  a 
strong,  static,  spatially  periodic  electric  Held  equal  to 
ft  Re{£  exp(ik  •  x)J.  This  field  in  turn  modulates  the 
index  of  refraction  by  the  first-order  electro-optic 
( Pockels)  effect  to  create  a  refractive-index  grating  in 


the  crystal.3  A  third  reading  beam,  also  at  u,  having  a 
wave  vector  k3  ■  — ki,  scatters  from  this  grating  to 
create  a  fourth  signal  beam  of  wave  vector  k+  »  —  k2, 
which  is  a  phase  conjugate  of  the  second  beam.4-5  (See 
Fig.  1.) 

Let  I\,  I %  1 3,  and  /4  be  the  incident  intensities  of  the 
writing  reference  beam,  the  writing  image  beam,  the 
reading  beam,  and  the  output  intensity  of  the  phase- 
conjugate  signal  beam,  respectively.  Consider  the  case 
in  which  all  four  beams  are  confined  to  the  y-z  plane, 
with  the  z  direction  taken  along  the  c  axis  of  the  crystaL 
(See  Fig.  1.)  According  to  our  previous  theory*  of 
grating  formation  in  BaTiOa,  when  /4  tS  0.5/;;,  the  mir¬ 
ror  reflectivity,  here  defined  as  ratio  R  »  IJI-i  of  the 
intensities  of  the  phase-conjugate  beam  to  the  image 
beam,  is  well  approximated  by 


£ord 
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n„3rl3  cos# 
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for  a  reading  beam  with  ordinary  polarization  and  by 
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Fig.  1.  The  reference  writing  beam  ( 1 )  and  the  image  writing 
beam  (2)  interfere  in  a  crystal  of  undoped  BaTiO;i  to  make  a 
refractive-index  grating  with  wave  vector  k.  The  reading 
beam  (3)  Bragg  scatters  off  this  grating  to  produce  the 
phase-conjugate  signal  beam  (4).  The  crystal  is  immersed  in 
index-matching  oil. 
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for  a  reading  beam  with  extraordinary  polarization  in 
the  crystal.  These  expressions  do  not  include  the  ef¬ 
fects  of  two- beam  energy  coupling6*6  or  the  phase  mis¬ 
match  that  are  due  to  the  change  in  the  index  of  re¬ 
fraction  from  each  beam  alone.*'  Here,  c  is  the  speed 
of  light  in  vacuum,  rj  a  (I.\H>Vri,  n„  and  nr  are  the  or¬ 
dinary  and  extraordinary  indices  of  refraction  in  the 
crystal  ( n„  =»  2.488  and  n ,  =  2.424  at  514  nm),  and  n.t  is 
the  index  of  refraction  for  the  reading  beam.  L  is  the 
effective  interaction  length  and  is  approximately  l 
expt  -  Vvyl),  where  l  is  the  beam  length  in  the  crystal  and 
7  is  the  optical  intensity-attenuation  coefficient.  The 
r„  are  the  conventional  contracted  electro-optic  coef¬ 
ficients  and  in  BaTiOj  (in  units  of  10-12  mV)  (Ref.  10) 
are  rt:l  =  8,  rXi  =  23,  and  r42  =*  820.  In  Eqs.  (1)  and  (2), 

H  is  the  angle  between  the  grating  wave  vector  k  and  the 
direction  of  the  c  axis,  and  aj  and  ct2  are  the  angles 
formed  by  each  writing  beam  with  the  y  axis.  From 
Ref.  6  one  sees  that  the  intensity  dependence  of  Eij  is 
contained  in  a  factor  fi1/2/s1/2/(f i  +  /•>  +  h)  so  that,  for 
I>  « 1 1  or/) «  7:  i,  the  mirror  reflectivity  ft  given  by  ei¬ 
ther  Eq.  (1)  or  Eq.  (2)  is  independent  of  the  incident 
intensity  /•>  and  depends  only  on  the  relative  intensity 
I:\JI  i  of  the  counterprsrpagating  beams.  Otherwise,  the 
electric-field  amplitude  E  depends  only  on  the  tem¬ 
perature  of  the  crystal  lattice,  the  charge  and  density 
of  the  migrating  carriers,  the  dc  dielectric  constants  of 
the  crystal,  and  the  relative  orientation  of  the  crystal 
and  the  optical  beams.6  In  general,  the  writing  image 
beam  will  form  an  intensity-interference  pattern  not 
only  with  the  reference  beam  but  also  with  the  reading 
beam.6  For  the  range  of  beam  angles  used  below,  we 
estimate  that  this  grating  contributes  about  10%  to  the 
observed  reflectivity  R. 

Inspection  of  Eq.  (2)  shows  that,  for  a  large  range  of 
angles,  the  reflectivity  for  extraordinary  beams  can  be 
larger  than  unity,  owing  to  the  contribution  from  the 
unusually  large  r42  coefficient.  (Previous  experiments 
with  BaTiOn  had  k  parallel  to  the  crystal  c  axis,  making 
t)  =*  0  and  thereby  precluding  any  contribution  from  the 
r42  term.)  For  example,  with  L  =s  0.4  cm  (7 1  «  1),  7j 
=  /  .  h  *  22°,  «i  =  18°,  «•>  -  26°,  and  a  calculated11 
value  of  E  *  4.4  X  102  V/cm,  we  compute  a  mirror  re¬ 
flectivity  of  ft ext  =  3.2  for  extraordinary  polarizations 
at  514.5  nm.  By  approximating  these  conditions  in  an 
experiment  with  I>  *  0.3  mW  and  Ii  =  /;t  =  5  mW 
(beam  area  —0.25  mm2),  we  observed  ftMt  —  2.  How¬ 
ever.  in  this  instance  we  used  ordinary  rays  for  beams 
1  and  2  to  write  the  grating  and  used  an  extraordinary 
ray  (beam  3)  to  read  the  grating  [which  does  not  alter 
Eq.  (2)  if  ki  —  k2  *  lt»  -  k;4}.  When  we  used  extraor¬ 
dinary  rays  for  all  three  incident  beams,  we  observed 
even  higher  reflectivities  (ft,xt  —  100)  because  of  the 
added  contribution  of  energy  coupling  between  beams 
1  and  2.  (With  the  above  geometry,  this  coupling  is  less 
than  0.02  for  ordinary  rays  but  can  exceed  50  for  ex¬ 
traordinary  rays,  and  the  coupling  greatly  enhances  the 
reflectivity  by  increasing  the  intensity  of  the  image 
writing  beam  as  it  propagates  through  the  crystal.) 

The  optical  setup  that  we  used  to  obtain  these  large 
reflectivities  is  shown  in  Fig.  1.  The  optical  beams  are 
incident  upon  the  barium  titanate  at  a  glancing  angle 
to  the  surface  of  the  crystal.  The  crystal  is  immersed 
in  index-matching  oil  (n  -  1.51)  in  order  to  increase  the 


Fig.  2.  A  plot  of  the  measured  mirror  reflectivity  /?,„  =  l  JI ■> 
as  a  function  of  the  reading-beam  intensity  I...  The  object- 
beam  intensity  was  fixed  at  /■_.  =  / 1/4,  and  the  angles  of  the 
incident  beams  were  «i  =  16°,  «•>  =  24°,  and  D  »  20®  (see  Fig. 
1 ).  In  this  plot,  I.\  has  been  normalized  by  the  fixed  intensity 
I-,. 


Fig.  3.  Optical  setup  for  observing  cw  self-oscillation.  The 
incident  beams  1  and  3  are  both  linearly  polarized  in  the  plane 
of  the  figure  and  are  extraordinary  rays  in  the  crystaL  Self¬ 
oscillation  is  observed  to  grow  between  the  crystal  and  a  94% 
reflectivity  plane  mirror  M.  Here  L’s  are  lenses  with  a  focal 
length  F  =*  100  mm.  and  P  is  a  variable  pinhole  used  to  control 
the  transverse-mode  structure  of  the  oscillation.  The  phase 
aberrator  A  is  formed  by  bubbles  of  transparent  glue  on  a 
microscope  slide.  The  angles  of  the  beams  are  about  the  same 
as  in  Fig.  2. 

angle  between  the  c  axis  and  k.  The  argon  laser  pro¬ 
duces  a  TEM,h>  Gaussian  mode  at  514.5  mm  in  a  single 
longitudinal  mode. 

The  following  four  experiments  elucidate  the  mirror 
characteristics.  In  the  first  experiment,  the  intensities 
/ 1  and  /2  of  the  writing  beams  are  fixed,  and  the  scat¬ 
tering  efficiency  of  the  grating  is  measured  as  the  in¬ 
tensity  I. s  of  the  reading  beam  varies.  All  the  beams  are 
extraordinary  rays  in  the  crystal.  From  Fig.  2,  it  is  seen 
that  the  intensity  of  the  reflected  beam  (or  signal  beam) 
can  be  made  to  exceed  the  intensity  of  either  of  the 
writing  beams.  Since  this  signal  beam  is  the  phase 
conjugate  of  one  of  the  writing  beams  (beam  2  in  Fig.  1). 
the  system  acts  as  a  phase-conjugate  mirror  with 
gain. 

In  the  second  experiment,  two  extraordinary  coun- 
terpropagating  beams  (beams  1  and  3)  are  incident  upon 
the  crystal  (beam  2  is  blocked).  A  plane  mirror  is 
placed  within  view  of  the  crystal,  with  the  normal  to  the 
mirror  directed  approximately  toward  the  crystal.  T wo 
new  counterpropagating  phase-conjugate  beams  are 
observed  to  grow  between  the  crystal  and  the  mirror 
with  a  time  constant  of  the  order  of  1  sec.  If  the  mirror 
is  tilted,  these  new  beams  will  fade  away,  only  to  reap¬ 
pear  on  whichever  part  of  the  mirror  is  closest  to  normal 
to  the  crystal.  In  the  process  of  finding  the  cavity  mode 
with  the  least  loss  in  which  to  oscillate,  the  crystal  finds 
and  directs  a  beam  at  the  most-reflective  surface  facing 
it.  Oscillator  output  fades  slowly  (  —  1  sec)  if  beam  1  is 
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blocked  but  extinguishes  instantly  if  beam  3  is  blocked, 
since  beam  1  is  helping  to  write  the  grating  but  beam  3 
is  reading  it. 

In  the  third  experiment,  two  identical  lenses  (L)  and 
an  aberrator  (A)  are  placed  in  the  cavity  formed  by  the 
phase-conjugating  mirror  (i.e.,  the  crystal)  and  the 
mirror  M,  as  in  Fig.  3.  Self-oscillation  is  allowed  to 
build  up,  and  its  transverse-mode  structure  is  photo¬ 
graphed  near  both  the  real  mirror  and  the  wave-front- 
reversing  mirror.  Figure  4  shows  these  mode  patterns, 
both  of  which  are  severely  distorted  by  the  aberrator. 
When  a  pinhole  is  placed  at  the  focal  length  of  the  lens, 
in  the  manner  suggested  by  AuYeung  et  al.,u  the  mode 
pattern  becomes  uniform.  If  the  crystal  were  acting  just 
as  an  ordinary  mirror,  the  light  distribution,  returning 
to  the  pinhole  from  this  mirror,  would  be  doubly  dis¬ 
torted  from  having  passed  through  the  aberrator  twice 
and  would  spill  out  and  be  blocked  by  the  face  of  the 
pinhole,  causing  a  large  loss  in  the  resonator.  In  fact, 
we  observed  that  when  the  adjustable  pinhole  is  made 
sufficiently  small  (800-um  diameter)  so  as  to  reject 
high-order  modes,  the  intracavity  power  in  the  reso¬ 
nator  decreases  at  most  by  about  10%  and  sometimes 
increases,  indicating  that  little  light  is  lost  on  the  walls 
of  the  pinhole  and  that  the  crystal  is  acting  as  a  high- 
quality  phase  conjugator. 

In  the  fourth  experiment,  an  image  amplifier  with  an 
intensity  gain  of  —10  is  constructed  by  using  two 
counterpropagating  beams  and  an  object  beam,  all  ex¬ 
traordinary  rays.  When  a  resolution  chart  is  placed  in 
the  object  beam,  an  amplified  real  image  of  the  resolu¬ 
tion  chart  is  observed.  (See  Fig.  5.) 


(»)  (b) 


(c)  (d) 


Fig.  4.  Photographs  of  far-field  mode  patterns  of  self-os¬ 
cillation  with  a  severe  phase  aberrator  in  the  resonator  cavity. 
(See  Fig.,'!.)  With  no  aperture  in  the  resonator  cavity:  mode 
pattern  transmitted  (a)  through  the  hack  of  the  crystal  and 
(hlthrough  the  94rt  mirror.  With  a  1-mm-diameter  pinhole 
in  the  cavity:  mode  pattern  (c)  from  the  back  side  of  the 
crystal  and  <d>  transmitted  through  the  94^  mirror.  These 
mode  pattern-  were  displayed  on  a  white  card  2  m  from  the 
cavity,  photographed  with  Kodak  Plus-X  (ASA  1251  film,  and 
printed  on  high-contrast  patter. 


Fig.  5.  Photograph  of  the  real  image  of  a  resolution  lest  chart 
(wheel  diameter,  1  cm)  formed  in  the  object  plane.  The  in¬ 
tensity  of  the  image  beam  was  measured  to  be  —10  times  the 
intensity  of  the  object  beam,  demonstrating  optical  image 
amplification.  The  bright  spot  seen  on  the  left-hand  side  is 
from  self-oscillation  between  the  crystal  and  one  of  the  faces 
of  the  glass  cuvette  that  holds  the  crystal.  This  image  was 
photographed  with  Kodak  Plus-X  (ASA  125)  film  and  printed 
on  high-contrast  paper.  The  angles  of  the  beams  are  about 
the  same  as  in  Figs.  2  and  3. 


In  conclusion,  we  have  demonstrated  a  cw  phase- 
conjugating  mirror  with  gain  up  to  100.  We  have  used 
this  mirror  to  construct  an  image  amplifier  and  an  op¬ 
tical  resonator  that  self-oscillates.  With  the  aid  of  a 
spatial  filter,  this  oscillating  resonator  will  correct  phase 
aberrations  inside  the  resonator  cavity  and  emit  a 
TEMih)  Gaussian  mode. 
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In  Fig.  1  the  angles  and  should  be  measured 
from  the  z-axis  instead  of  from  the  y-axis  as  shown. 
Therefore,  the  17th  and  18th  line  on  page  520  should  be 
changed  to  read  "....  and  and  c*2  are  the  angles  formed 
by  each  writing  beam  with  the  z  axis.  From  ..."  On 
the  42nd  line  of  the  same  page,  the  values  given  for 
and  o2  should  each  be  increased  by  90°,  as  should  the 
values  of  and  a2  given  in  the  caption  of  Fig.  2.  This 
correction  does  not  alter  any  of  the  calculations  or 
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